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Gaussian and non-Gaussian fluctuations for mesoscopic linear 
statistics in determinantal processes 

Kurt Johansson* Gaultier Lambert^ 


Abstract 

We study mesoscopic linear statistics for a class of determinantal point processes which in¬ 
terpolates between Poisson and Gaussian Unitary Ensemble (GUE) statistics. These processes 
are obtained by modifying the spectrum of the correlation kernel of the GUE eigenvalue process. 
An example of such a system comes from considering the distribution of non-colliding Brownian 
motions in a cylindrical geometry, or a grand canonical ensemble of free fermions in a quadratic 
well at positive temperature. When the scale of the modification of the spectrum of the GUE 
kernel, related to the size of the cylinder or the temperature, is different from the scale in the 
mesoscopic linear statistic, we get a central limit theorem (GET) of either Poisson or GUE type. 
On the other hand, in the critical regime where the scales are the same, we get a non-Gaussian 
process in the limit. Its distribution is characterized by explicit but complicated formulae for the 
cumulants of smooth linear statistics. These results rely on an asymptotic sine-kernel approxima¬ 
tion of the GUE kernel which is valid at all mesoscopic scales, and a generalization of cumulant 
computations of Soshnikov for the sine process. Analogous determinantal processes on the circle 
are also considered with similar results. 
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1 Introduction and results 

1.1 Introduction 

Recently there has been a lot of discussion about universality of random matrices statistics at mesoscopic or 
intermediate scales. For instance, the proofs of the local semicircle law and the Wigner-Dyson-Gaudin-Metha 
conjecture, see [mill US] for further references, the work [miiH] on random band matrices and the so-called 
Anderson metal-insulator transition, or the GET for linear statistics of orthogonal polynomial ensembles |8] . 
One motivation to investigate these models comes from E. Wigner’s fundamental observation that the spectral 
statistics of complicated quantum systems exhibit universal patterns. On the other hand, eigenvalues of 
quantum systems which are classically integrable are expected to be described by Poisson statistic m- 
Therefore, it is natural to investigate the transition from Poisson to random matrix statistics at intermediate 
scales. There are many ways to interpolate between the two point processes, for instance using Dyson’s 
Brownian motion, one gets a determinantal process called the deformed Gaussian Unitary Ensemble |24j . For 
this model, the transition has been investigated using mesoscopic linear statistics in M- The authors proved 
some GLTs with fluctuations which depends on the scale the test function samples the eigenvalues. In this 
paper we will study the fluctuations of another general class of determinantal processes which interpolate 
between Poisson and GUE statistics that we call modified GUEs (see definition II.ip . Instead of adding an 
independent matrix to the GUE, we directly modify or mollify the spectrum of the correlation kernel of the 
process. This has the effect of introducing some extra disorder in the system while keeping the determinantal 
structure. Our main motivation to study such ensembles comes from the so-called MNS ensemble which was 
introduced in |28l[24] and corresponds to the joint p.d.f. of the positions of a Grand-Canonical system of Free 
fermions at positive temperature confined in a one dimensional harmonic potential [10] . In general, it is also 
of interest to investigate fluctuations of determinantal processes whose correlation kernels are not necessarily 
reproducing. 


1.2 The Modified Gaussian Unitary Ensembles 

Let X be a Polish space equipped with some reference measure dfi. In the sequel, we will only be interested in 
the two cases X = M or the unit circle T equipped with the Lebesgue measure. A point process is as a random 
measure on X of the form S = ^ 6xi- The support of the measure S is the random object of interests, it is 
called a point configuration (Xi) and we assume it has neither double point nor accumulation point. Point 
processes are usually described by their correlation functions pk{xi,... ,Xk) which are characterized by 


E 


n(i+5(^*)) 


[ U9{x^)pk{xu 


.. .,Xk)dp,{xi) ■ ■■dp,{xk) 


( 1 . 1 ) 
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for any measurable continuous function g : X —>■ C with compact support. A point process is called determi- 
nantal if its correlation functions exist and satisfy the identity 

Pkixi, ...,Xk) = det[K{xi,Xj)] . 
kxk 

Hence a determinantal process is characterized by its correlation kernel K : X x X —>■ M and we will denote 
by E/f the corresponding probability measure on the space of point configurations. One generally assumes 
that K defines an integral operator K, on (X, dp) which is locally of trace class and in such case equation 
dm can be seen as the Fredholm determinant, 


Ek 


n(l +5(a:*)) 


— det[/ + . 


( 1 . 2 ) 


In most cases the operator K, is self-adjoint (although there are natural examples of non-Hermitian deter¬ 
minant processes, such as the deformed GUE studied in [2]). Then, the kernel K defines a determinantal 
process if and only if all the eigenvalues of the operator K. lies in [0,1]. These facts are well-known and we 
refer to [221E31131] for different introductions to the theory of (determinantal) point processes and to the 
survey [3] for an overview of some applications. In this paper we will investigate examples of determinantal 
processes with correlation kernels of the general form 


K{x,y) = ^p^'dkix)'dkiy) , 
A;=0 


(1.3) 


where (dfe) is an orthonormal basis in L^{X, dp) and the spectrum 0 < p^ < 1. 

A classical example of such correlation kernel is that of the GUE eigenvalue process 

N-l 

Kq {x, y) = Y^ ipk{x)ipk{y) , (1-4) 

k=0 


where pkix) = ^.re the rescaled Hermite functions, see (15.11) below. We refer to [53] 

Ghap. 2-5 for some background on the GUE process. Let us just mention that the parameter N corresponds 
to the number of eigenvalues and the process is scaled so that its density at the origin is equivalent to N as 
N ^ oo. We are mainly interested in determinantal processes whose correlation kernels are modihcations of 
. Namely instead of taking the spectrum p^ = lk<N, we assume that k i—>■ p^ is a function which decays 
from 1 to 0. 


Definition 1.1. A modified GUE is a determinantal process on K (wrt. the Lehesgue measure) whose 
correlation kernel ^ depends on a function £ 5^ called the shape, see (11.61) . and three parameters 
N > 0, a € {0, 1) and t > 0 in the following way 

K^,aix,y) Ma ] ^k{x)pk{y) , (1-5) 

where (fikix) = \Jhk Hermite functions given by (15.11) . In the sequel, the 

parameter a is called the modification scale and t the temperature. 

We will consider the following class of shapes, 

S' = {'h : R —>• [0,1] I T' < 0 is Riemann integrable, 'k £ L^(0, oo), (1 — '!') £ L^(—oo, 0)} . (1.6) 

If T £ S) one can think of 1 — 'h as the distribution function of some probability measure on R and of 
$ = —'L' as the corresponding density. For technical reasons, we will also consider the following subclass 

S' = £ S I 4>(x) = —'I''(a:) < for some c > 0} . (1.7) 
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The kernel ^ is not reproducing, so the total number of particles in the process, denoted by is random. 
Moreover, the process has been scaled so that its density at the origin is ^ N (global scaling) and a simple 
computation yields ^ N. According to theorem 7 in [22], an equivalent correlation kernel for such 

process is given by 

OO 

K{x,y) , (1-8) 

k=0 

where are independent Bernouilli random variables with E [l^] = ^ ) ■ Hence the modified en¬ 

sembles are more random than the GUE and in some sense the amount of extra randomness is estimated 
by 

OO 

Var^«J#] = ^Var[/f] . (1.9) 

A Riemann sum approximation gives 

Var^iv [#] - rAf“ [ ^'(t)(l - ^-(1)) dt . (1.10) 

Jr 


Heuristically, it means that the more the spectrum of the correlation kernel is modified, the more disorder is 
forced into the system. So we expect that, at large modification scales, the modified ensembles behave like 
the Poisson process rather than like the GUE. Note that, by (II.IOL the assumption G 5" guarantees that 
Var^fW [4f\ < OO for any iV > 0. 

Our interest in determinantal processes with correlation kernels (HU is mainly motivated by the following 
example that we call the MNS ensemble. In the physics literature, the authors of [55] introduced an 
ensemble of unitary invariant Hermitian matrices whose eigenvalue distribution interpolates between the 
GUE and the Poisson process. This model was rigorously analyzed in [53] and it was proved that its Grand 
Canonical version is a determinantal process with correlation kernel on R given by 


irN f 


y) = 


1 


1 + eik-N)/rN- 

k—0 


'(Pk{x)(pk{y) 


( 1 . 11 ) 


for some r > 0 and 0 < a < 1. So the MNS ensemble is a modified GUE with shape ipit) = (1 -I- 
This model has two natural interpretations. First, since the Hermite functions are eigenfunctions of the 
Schrodinger operator — A -|- on R., the process describes a grand canonical system of free fermions at 
positive temperature confined in a quadratic external potential. Namely, the probability that the k**^ state 
of this harmonic oscillator is occupied is equal to the Fermi factor (1 -|- where (3 = (A“r)“^ is 

the inverse temperature of the system. The Gaussian Unitary Ensemble corresponds to the ground state of 
such a system with N fermions. Namely, taking the temperature to zero (i.e. the limit r —>■ 0 in (11.111) 1. one 
recovers the GUE kernel given by (HH). On the other hand, for large temperature (i.e. taking r —>■ oo), 
the kernel degenerates to that of a Poisson process on R. Therefore, at a heuristic level, the MNS ensemble 
interpolates between Poisson and random matrix (GUE) statistics. We shall prove that such a transition 
occurs for smooth mesoscopic linear statistics of the process. Second, in [24], it was shown that the MNS 
process also describes a system of Brownian particles moving on a cylinder and conditioned to never intersect 
(by rotation invariance, the distribution of the particles is stationary). We have seen that the parameter N is 
the expected number of particles and one can check that in this case the parameter P = corresponds 

to the perimeter of the cylinder around which the particles are diffusing. This process behaves roughly like 
Dyson Brownian motion UM and this provides a heuristic description of the transition. Namely, at small 
scales, the particles remain roughly independent, while when /3 gets large the trajectories start regularizing 
because of the non-intersecting constraints until eventually their joint distribution obeys the law of the GUE 
eigenvalues. 
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In random matrix theory, it is well-known that we can analyze the eigenvalues processes at different scales. 
The global or macroscopic scale refers to the size of the whole process. On the other hand, the local or 
microscopic scale is that of individual eigenvalues, i.e. the gap between two eigenvalues is of order one. At 
this ultimate scale, in the Hermitian case, universality roughly means that the process converges in the bulk 
to the celebrated sine process. Any scale in between is called mesoscopic. In other words, a mesoscopic 
random variable is a function of the point process which depends on a growing fraction of the total number 
of particles. A typical example of such observable is that of mesoscopic linear statistic. 

Definition 1.2. Given a point process S with density ^ N at the origin and a function / : X —>■ C with 
compact support, for any 0 < 5 < I, we define 

fs{x) = f{xN^) , 

and we call a mesoscopic linear statistic the random variable 

Efs = Y,f{X,N^) . ( 1 . 12 ) 


In the sequel, the parameter 6 is called the scale. 

With such convention, the scale i5 = 1 is microscopic, the scale 5 = 0 is macroscopic, and we write S/ for 
S/o. In the sequel, we will investigate the distribution of mesoscopic linear statistics of the modified GUEs. 
Since the density at the origin is of order N, 


[Efs] ^ [ f{x)dx . 

’ Jr 


(1.13) 


If <5 < 1, this expectation is diverging as N ^ oo and it is natural to consider centered linear statistics instead 

Efs = Y. [S/^] • (1-14) 

i 

For a determinantal process, there is an explicit way to express the cumulants of the random variable 
S/ in terms of the correlation kernel AT; see formulaensemble as the density N ^ oo. 


1.3 Main results 

Throughout this paper we will use the notations from definitions 11.11 and 11.21 the parameter S G [0,1] will 
always refer to the scale of a linear statistic and the parameter a G (0,1) to the modification scale of the 
correlation kernel (USD. All other formulae and definitions are collected in section 11.61 In theorems 11.4111.61 
we assume that the shape of the correlation kernel has an exponential tail, i.e. 'k G see (11.71) . This 
assumption is not necessary but it does prevent some technical difficulties; see remark [SA] and the discussion 
after condition (14.511 . Moreover, it applies to the MNS ensemble which is our main example. Our results are 
summarized in the diagram of figure [TJ 

Theorem 1.3. For any parameters 0 < 5 < a < 1 and for any function f G with compact support, 

Var^. [SM = ^iV“-^B| //(x)2dx+ o (fV“-^) , (1.15) 

where the shape-dependent constant 0 < < oo is given by (11.411) . This asymptotic formula implies the 

following classical Central Limit Theorem as N ^ oo, 

N-^ Efs ^ Ar(0,2rB|||/||i.) . (1.16) 
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Proof. The asymptotic expansion of the variance is proved in section [ST] The CLT (11.1611 follows directly 
from the estimates (11.131) and (11.151) by applying Theorem 1 in [3B]. □ 


Hence, we will call the set S [0,1) : 5 < a] the Poisson scales because the variance of linear statistics is 
diverging in this regime. Viewing the point process 2 as a random distribution in iS'(]R), theorem II.31 implies 
that a rescaled modified GUE converges at any scale <5 < a to the white noise with intensity tB^, i.e. a 
Gaussian field 2oo on the real line with covariance structure 

E [2oo(/)Soo(ff)] = TB|(/,5)i2(R) . 

The characterization of the process at scales 6 > a is more subtle because the fluctuations remains bounded. 
Therefore we expect a limiting process with strong correlations. We will prove that linear statistics of the 
modihed GUEs are close to that of the sine process at any scale S > a, and from this deduce the following 
GET. 

Theorem 1.4. Let 4' G 5^' and f € with compact support. For any parameters 0 < a < 5 < 1, as 

N —>■ oo, 

Efs ^ Af (0, 11 / 11 ^ 1 / 2 ) , 

where the norm ||/||^i /2 is given by formula (11.281) . 

Proof. Section 13.21 □ 


Hence, we will call the set G (0,1) : a < 5} the GUE scales by analogy with theorem l3.6l An interpretation 
of theorem ll.ll is that the centered modified GUEs converge weakly at any scale i5 > a to a Gaussian process 
2o on the real line with covariance 

E [2o(/)2o(5)] = {f,g) hi /2 . 

Gontrary to the white noise, the Gaussian process 2o is spatially correlated and self-similar as can be seen 
from equation (11.281) . 


Theorems o and HH imply that any modified GUE undergoes a transition from Poisson to random matrix 
statistics when the mesoscopic scale S is equal to the modification scale a of its correlation kernel. Our next 
question is what happens at the critical scale? The first step is to compute the limit of the variance of linear 
statistics. Depending on the context, we obtain different formulae for the critical variance and we check 
that they are consistent in appendix [^ Theorem 11.51 summarizes our results and it is also a special case of 
theorem 11.61 

Theorem 1.5. For any shape T G , any function f G with compact support, and any scale 

0 < a <1, 


lim Var^N [^fa] = 2 t' 

N^oo 


= 2t' B| 


' f f{u) f T'(t)(l — T'(t-I-u/r')) dudt 

Js. JS. 


fix) - fiy) 


x-y 


$(r'(x - y)) 


dxdy 


(1.17) 


where $ = — tk' and the parameter P = j > 0. 

Proof. Appendix]^ □ 


Since $ is a p.d.f., it is clear that formula (11.171) interpolates between ||/||^i /2 as t —>• 0, respectively 
2tP\ II/III 2 as r — 00 . In both cases, we recover the variance of theorem 11.31 respectively theorem 11.41 We 
would expect Gaussian fluctuations when ^ = a, in analogy with the other cases. Surprisingly the cumulants 
of linear statistics of the modified ensembles have non-trivial limits at the critical scale and thus we get 
a non-Gaussian limit. In order to formulate our main result, we need to introduce several notations; see 
definitions (11.381) and (11.391) in section [T^ 
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Theorem 1.6. Let f € with compact support, 0 < a < 1 and G . The linear statistic Sfa 

of the determinantal process with correlation kernel given by (11.51) . converges in distribution as the 

density N ^ oo to a random variable ^^ 4 /. Its cumulants are given by 


j fUTdt -2 Y, M{-m) JJ ^lf[f{u,)^xAG^{u,x)d^-^ud^x . (1.18) 

® |m|=ra R^xRI ' 


Proof. Section 14(2 


□ 



Figure 1: Fluctuations of a Modified GUE as a function of the scale <5 and the modification length a. 


Unfortunately, it seems very difficult to recover from equation (11.1811 the Laplace transform of or to 

say much about its distribution. However, we can prove that, provided there exists n > 2 such that 0, 

the random process is not Gaussian; see proposition 14.141 and (11.391) . This condition is very general, 
but somewhat surprisingly it does not apply to the MNS ensemble. In fact we have the following result. 

Proposition 1.7. The shape function ^^{t) = (1 + satisfies = 62 {n). 

Proof. Section 14.31 □ 


In particular, by equation (I4.65L the conclusion of proposition 1 1. 71 is equivalent to the fact that the function 
Ip solves the equation 


/ 


-dt = 


VfeC 


(1.19) 


1 + 'ip{t)(e^ — 1) e? — 1 

This shows that the MNS ensemble is a special process among the class of modified GUEs. This property 
is maybe related to the fact that this point process originates from a Grand Canonical model and it would 
be interesting to know whether equation (I1.19P has any other solution in The main consequence of 
proposition ll.7l is that we need a separate argument to prove that the limit of the critical MNS ensemble 
is not Gaussian. To this end, we show in section im that the function y{x) = cos(27ra;) G 5(]R) satisfies 

G" [E^,ry] 0 , (1.20) 
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when e is sufficiently small. We need to look at the 4**^ cumulint since = 0 for any test function 

because of some symmetries; see proposition 14.191 There is nothing special about the function y{x) except 
that it is simple enough to provide a good example. Equation (11.201) should hold for most test function but 
this is certainly very technical to prove. 

1.4 Discussion and Examples 

We have seen that the MNS ensemble describes non-colliding Brownian motions in a cylindrical geometry. In 
this case, the diagram of figure [T] shows the particles statistic at a given mesoscopic scale 5 exhibits a sharp 
transition from Poisson to GUE at time This situation is similar to Dyson’s Brownian motion which 

was investigated in [T3]. However the transitions are different: for the MNS ensemble there is no interme¬ 
diate regime which depends on the test function / and unexpectedly the fluctuations are not Gaussian at 
the critical scale. In both interpretations of the MNS model as non-colliding Brownian motions, or a grand 
canonical ensemble of free fermions, it is not clear why we get non-Gaussian fluctuations at a certain relation 
between the sampling scale of the linear statistic and the size of cylinder, respectively the temperature of 
the gas. Actually, it would be very interesting to get another description than theorem 11.61 of the random 
field which arises at the critical scale in order to understand how the spatial correlations disappear in 
the transition from the iJ^^^-Gaussian noise Sq to the white noise Sqo- Another point of interest would 
be to understand how the processes depend on the shape function 4 and if proposition 11.71 has any 
probabilistic or physical interpretation. 

For determinantal processes, a general strategy to get a GET for linear statistics is to use equation (11.311) 
and show that all cumulants of order > 3 converge to zero. This approach has been applied to many classical 
models in 1 or 2 dimensions using a wide range of techniques, see e.g. [glEIllMlEIlEollIllTlIS]. In this paper, 
it is applied to the modified ensembles of definition ll.il To the authors’ knowledge, the results of this paper 
at the critical scale provides the first example of a determinantal point process for which the fluctuations of 
a mesoscopic linear statistic are not Gaussian in the limit. At the global scale non-Gaussian limits can be 
obtained for Hermitian matrix models with several cuts, but the mechanism there seems to be very different. 
To prove theorems 11.41 and 11.61 we use a perturbative approach which consists in comparing the correlation 
kernels of two processes to show that a given linear statistic S/ has the same limit for the two processes as 
the density N ^ oo. To this end we will use the following definition. Let be the n*** cumulant for 

the linear statistic S/ of the determinantal process with kernel K. 

Definition 1.8. Two families of kernels (Ar^)o<x and {L^)o<n defined on the same space 3C are asymp¬ 
totically equivalent {we write = L^) if for any function f £ C'o(X) and for all n G N, 

hm = hm Cl^Ef] . 

N—¥(x> N—^oo 

Note that definition II.81 does not only apply to correlation kernels since we only need to assume that a kernel 
K is locally trace-class to define C^lEf] according to equation p.3ip . However, if both and define 
determinantal processes and , then these processes have the same limit as A^ —>■ oo. For instance, 

theorem 11.41 is proved by showing that the kernel ^ and the GUE kernel Kq , see (11.51) resp. (11.41) . are 
asymptotically equivalent at any scale 6 > a and using the GLT for the mesoscopic GUE (theorem [32] proved 
in [SlISlEnilllEe]). On the other hand, at the critical scale (i5 = a), the kernels ^ are not asymptotically 
equivalent to any kernel which has been studied previously and we will need to compute the limits of the 
cumulants explicitly. If 4 G S’, we let and define G ( 0 , 1 ) such that by the mean-value theorem 

7^^ ~ (t^) ~ fc G Z. Then, if A^, r, T > 0, a G (0,1), and 77 is a non-decreasing 

function, we define the kernel 


y) 


1 

T-Ar« 


E 

|fe|<rAr' 


$ 


+ Cfe 


sin [27r?7(fc)(a: — y)] 


tN° 


n{x - y) 


( 1 . 21 ) 














Lemma 14.11 shows that in general the kernel ^ defines a determinantal process on K. and the following 
proposition implies that, for a given function t], linear statistics of the critical modified GUE with shape 4/ 
and the process with kernel ^ have the same limits. 

Proposition 1.9. Let G 5^' and 1/3 < a < 1, if 'q{k) = “-^/l + k/N /2 and L = (log the correla¬ 

tion kernel of the modified GUE satisfies N~‘^K^^{x,y) in the sense of dehnition W M 

Proof. Section 14.11 □ 


The main technical challenge to prove proposition 11.91 is to get a uniform asymptotic expansion of the GUE 
kernel valid at mesoscopic scale, theorem 15.31 fsee also formula (I5.1ip l. The restriction 1/3 < a < 1 comes 
from the fact that at larger scales, the GUE kernel is not asymptotically translation-invariant because of the 
curvature of the density of the semicircular law. In this case, we can still use the formula of theorem 15.21 
and make changes of variables in order to use the kernel ^ to compute the limits of the cumulants of the 
critical modified GUEs; see proposition 14 .1 ll Since the kernel L^ ^ is translation-invariant, we can compute 
the cumulants of its linear statistics by using the method introduced in [33] to prove a GLT for mesoscopic 
linear statistics of the Circular Unitary Ensemble (see also theorem 4 in [35] for an application to the sine 
process, as well as some generalizations). When we pass to the limit as ^ oo, we get theorem 11.61 In 
fact, taking the parameter r —>■ 0 in equation (11.181) . we recover Soshnikov’s formula for the mesoscopic sine 
process: 


C" 





M(m)G™(u) d^-^u 


where G™(m) = ^ max{u^(l)H- \-u^(rni), ■ ■ ■, u^ii) + ■ ■ ■ + O} ■ 

7rGS(n) 

Then, the Main combinatorial Lemma of |3S] implies that ^^M(m)G™(M) = 0 for any n > 2 and that 
the process (which is independent of 4') is Gaussian, as we expected from theorem 11.41 The details 
are given in the proof of proposition 14.151 For the modified ensembles, there is no counterpart of the Main 
combinatorial Lemma, it means that in general M(m)G™(u, cc) ^ 0 for any r > 0 and n ^ 3. This fact 
is used to prove equation (ll.20[) ; see in particular lemma 14.211 This shows that the combinatorial structure 
behind the cumulants of the sine process, which corresponds to the Strong Szegd theorem, is very sensitive. 
In general, CLTs with bounded variance are due to some special correlation structures which are rather 
sensitive under perturbation such as some small modification of the correlation kernel. In the remainder 
of this section, we provide a simple example which elaborates on this fact and illustrate how asymptotic 
normality breaks down. Before proceeding, we define the Circular counterparts of the modihed GUEs. These 
point processes are of interest because asymptotic expansions are not required in order to apply Soshnikov’s 
method. In the sequel, we let T = [—i, i] with the boundary points identified. 

Definition 1.10. A modified CUE is a determinantal process on T {wrt the Lebesgue measure) whose 
correlation kernel is of the form 


feGZ 


( 1 . 22 ) 


where = p^ so that the corresponding integral operator is self-adjoint on L^(T). 

In the rest of this section and section \27i\ the spectrum of the kernel is arbitrary except for the constraint 
p^ G [0,1]. However, when we refer to a modified CUE afterwards, it is understood that 


Pk='^ 


\k\-N 

tN°‘ 


(1.23) 


for some shape dt G S’, a G (0,1) and t > 0 by analogy with definition W.W 
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(1.24) 


A special case is the so-called Dyson CUE which has the correlation kernel 


K^{x,y) = ^ e 

\k\<N 


i2nk(x-y) ^ sin((2Af-H)7r(a;-y)) 
sin(7r(a: — y)) 


This process describes the eigenvalues of a random matrix distributed according to the Haar measure on 
U{2N -I- 1). The cumulants of its linear statistics were computed explicitly in [35]. A simple of extension of 
this proof yields the following formula in the case of the modified CUEs. 


Lemma 1.11. For any continuous function / : T 


i)-i 


c^^[Ef]= Y[fM n Pk+u^+-+u^^, 

uGZq i |m|—n i—0 


where f(k) is the Fourier coefficient of the function f. 
Proof. Section \^fT\ 


□ 


We can use this expression to investigate the behavior of (global) linear statistics under some very simple 
modification of the spectrum of the CUE correlation kernel. Eor instance, we can remove a single mode, i.e. 
we let 

Pk = ^\k\<N — ^\k\=N-m (1-25) 

for some m G Note that in this case is still a projection kernel. By lemma IT. Ill the 3'''^ cumulant 
of a linear statistic of a modified CUE is given by 


/^3 r'=' 


/] — ^ ^ J_ J_ / (^i) > ^ ^ Pk n ^ '. Pk{pk+Ui F Pk+Ui+U2') F PkPk+uiPk+ui+U2 

feez 


uGZg i KkSl. fcGZ 

We can symmetrize this expression using permutations of the ufs and the condition ui F U 2 = —M 3 , we get 

^/]=En f (^g) ^ ^ Pk {f ^Pk+ui T 2p/,._|_^^p^_^2} ■ (l-2b) 


C^n[“ ' 


feez 


Consider the function gj{t) = 2cos(27rjt) -|- acos( 47 rjt) for some a G R and j G Z+, so that 

= Sj{\u\) F ^S 2 j{\u\) . 

With this test function, the only frequencies u G Zq which contribute to (11.261) are given by all possible 
permutations of (±j, ±j, =F2j). A simple computation shows that 




(1.27) 


fcGZ 


In the CUE case, when p^ = l|fe|<Ar, an easy computation gives Cf.N[E.gj] = 0 for any j G Z+.This was 
expected, however it seems clear that for some generic choice of coefficients 0 < p^ < 1, the expression (|1.27|) 
will be non-zero. Eor instance in the case (|1.25|) . it is easy to check that for any j, m N, 

C^N[^gj] = 12 a(l — lj<[m/ 2 j) ■ 

Hence, if we remove a mode near the edge of the spectrum of the CUE kernel (m < 2j), the linear statistics 
Spj is not Gaussian in the limit N ^ 00 . Moreover we can check that the variance is bounded with the same 
choices of parameters, 

CknI^Pj] = \gj{u)\^ YPki^ - Pk+u) < j{2Fa^) . 

u^Z k^Z 
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This implies that the limit as N ^ oo oi the determinantal process with correlation kernel is not a 
Gaussian process. In particular, even if the correlation kernel is reproducing, we can get non-Gaussian 
behavior. This example also shows that it is the edge of the spectrum of the correlation kernel which is 
influencing the distribution of the point process, see also theorem 4 in [36]. Moreover, by lemma IT. Ill we can 
check that when removing M modes at the edge of the spectrum, all cumulants are 0{M) for large N. Hence, 
if we remove sufficiently many modes, the system begins to behave like a Poisson process when iV —>■ oo. 
Finally, note that according to equation (|1.8I) . removing modes is comparable to smoothing the spectrum of 
the correlation kernel. Thus this example motivates why the modified ensembles of definition 11.11 are not 
Gaussian at the critical scale {5 = a). Actually, the strategy to prove equation (|1.20l) is the same as in this 
example but the computations are much more complicated. 


1.5 Overview of the rest of the paper 


In section |2l we begin by analyzing the modified GUEs of definition 11.101 This setting is simpler than that 
of section 11.21 and we can focus on the combinatorial structure of the cumulants. Lemma 11.111 is proved in 


section [2m and in section [2^2.31 we show that, if the spectrum of the kernel is given by = 4' 




then the results of figure [T] hold for the modified GUEs as well. The main results are proved in sections |3] 
and 01 The asymptotics of the variance in the Poisson regime, formula (11.15|) , is computed in section 13.11 
while in the GUE regime, theorem ll.dl is proved in sectionAt the critical scale, proposition [T!9| is proved 
in section im Both the critical modified GUEs and GUEs are analyzed, in a common framework, in section|4| 
In particular, in section [4.21 theorem 11.61 is proved in two steps. First the limits of the cumulants are given 
in theorem 14.41 (see also proposition 14.131 for large scales). Then the weak convergence of linear statistics is 
established in corollary 14.101 In section 031 we show that the random processes Sij, t- of theorem ll.61 are not 
Gaussian and we prove the special property of the MNS ensemble, theorem ll.71 by computing the generating 
function of the coefficient B^. We also prove that, as it is expected from figure lU the random field 
converges to a Gaussian process in both limits r —^ 0 and r —>• oo; see proDOsition l4.151 Finally, in section IU4l 
we show that the critical MNS ensemble is not Gaussian by giving an example. All the asymptotics that are 
required to analyze the modified GUEs are gathered in section |S| In appendix we prove theorem 11. 51 and, 
as an example, we compute the critical variance of linear statistics of the MNS ensemble. In appendixjB] we 
prove a technical lemma which is needed to obtain equation (11.20|) . Finally, all the notations that will be 
used in the rest of the paper are collected in section [TT6l 


1.6 Notations and cnmnlants formnlae 


We will denote 


and for any z € R, 


K" = {a; G M" : a;i < • • • < Xn} , 
R" = {x G K" : xi H-+ x„ = z} . 


The same definitions holds for Z", etc. 


We will also use the following conventions, 

• xat ^ zat if lim7v_>.oo xn /zn = 1 ■ 

• Xn — z_N if limAr_j.oo(xAr - z^) = 0. 

• Xn = 0 {zn) if there exists two non-negative constants k and C such that |x 7 v| < Czn\ logiVI”. 

Test fnnctions. 

For any function f G D L^(R), we define its Fourier transform 

f{u) = [ . 

Jr 
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We will consider the following spaces of test functions: 

C'o(X) denotes the space of continuous real-valued functions with compact support on X (any Polish space), 
denotes the Sobolev space of real-valued functions equipped with the norm 


Ih1/2 


[ fH 

Jr 


2 1 
\u\du=^ 


R 2 


fix) - fiy) 


x-y 


dxdy . 




denotes the Sobolev space of real-valued functions equipped with the norm 


\h^ 


Ku) 


l^du = 




dx . 


(1.28) 


(1.29) 


Modulo constants, the spaces i?^/^(]R) and H^{R) are normed spaces. We will also let Hq be the set of 
compactly supported functions in H'^. It is well known that C'q(K) C Hq(R) C i?g^^(M) C C'o(M) C 
n L^(R). Finally 5(R) is the space of Schwartz functions. 


Cumulants of linear statistics of determinantal processes 

For any random variable Z with a well-defined Laplace transform, its cumulants C'”[Z] are given by the 
formal power series 

j.n 

logE[e‘^] =^C"[Z]- . (1.30) 

Ti. 

n—1 

For determinantal processes, it turns out that there is a rather simple way to express the cumulants of a 
linear statistic S/ = fi^i) in terms of the correlation kernel. If we take g{x) = — 1 for some 

function / G C'o(X) in equation (11.21) . 


Ek 




det 


I + 



Since the operator K, is assumed to be locally of trace-class, the r.h.s. of this equation is well-defined and 
taking the logarithm, we obtain 


log Ek 




= Tr 


log (/ + /c(e‘^(^) -1; 


- (_iy+i 

“EM— Tr 


Z =1 



_^ -j-n 

(See for instance [33] chap. 3.) If we expand — 1 = 7 /(x)"— and use linearity of Tr, we get the 

n\ 

following expression for the cumulants of the random variable : 


C n 

K 


/i = E 


(- 1 ) 


Z +1 


E 


n\ 




mi >1 

mi-\ - \-7ni—n 


mi! ■ ■ ■ mil 


Tr[/™i/C/™" •• •/’"'/C] , 


(1.31) 


where we interpret as multiplication operators acting on L^(X, d/i). In the sequel, we will abuse notation 
and denote any integral operator /C like its kernel K. 


A composition of a number n G N is a tuple m = {mi, m 2 , ■ ■ ■ ,mi) of positive integers such that |m| = 
mi + • • • + m( = n, where i = f(m) is called the length of m. Denote the multinomial coefficient by 


n\ n\ 


and write 


M(m) 


(-ir+V nN 

i V™/ 


(1.32) 
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(1.33) 


By convention, we set (^) = 0 whenever |m| ^ n. Then we can rewrite equation (11.3111 as 

C^[Ef]= E M{m)Ti:[r^Kr^---r‘K] , 

\m\—n 

where 

Tr[r^Kr----r^K] = [ f{xl^^Kixl,X2)■■■f{xe^‘Kixi,xMxl)■■■d^i{xe). (1.34) 


Cumulants of the random variables 

For any composition m and for any k < i{ui), we define 

Wfe = mi + m2 + • • • + ruk ■ (1.35) 

The map m m can be viewed as a linear transformation on For any m S RI™!, we define 

rrli rus j ^rris + l 4“ * * * “t“ Ujn-i if S ^ i 

^is (^) ^ ^ ^ ^ \ + l ‘ ‘ if Z < S , (1.36) 

1=^ 4=1 1^0 if 4 = s 

so that A™ = (A™) is an ^(m) X ^(m) antisymmetric matrix. 


If CT G S(n) is a permutation of [n] = {1,..., n}, we will use the shorthand notation au = (ucr(i)) ■ • ■, Wcr(n)) 
and we define s{a) = argmin(cr) G [1] x [2] x • • • x [n] as follows. For any I — 1,... ,n, the number S/(cr) is 
given implicitly by the relation 

cr(s;) = min{cr(j) : j = 1,. . . J} . (1.37) 

For any r > 0 and any composition m of n > 2, let 

G^{u,x)= E max|A™ (u)-r(x^(q-a;<,(^^))| , (1.38) 

{7GS(n) ~ 

where we used the shorthand notations i = £(m) and Sg = S£(m)(o')- Note that, if cc G R", by definition 
> 2:^(8^) and A™j,^(u) = 0. Hence the functions G^{u,x) are non-negative on R” x R”. 


Finally, for any ih G 1?, let = 0 and for any n >2, 

n—1 


— 

-D^ — 


E^fe / a;$(x)«'(x)'=(l - «'(x))" ^ ''dx , 
k=0 •'* 


where $ = — ih' and the coefficients 6^ are given by 


fc-i-i 


^fc = E(-i)' 


i+i 


1=1 


n — I 
k + 1 — I 


E 

|m|=n 


(1.39) 


(1.40) 


These numbers are analyzed in section l473l For now, let us observe that = —b^ = 1 and by (11.391) . 

B^ = J a:$(a;)(l — 2il'(a;))da; 

Using that <i)(a:)(l — 2'l>(a:)) = —^{iE'(a;)(l — il'(a:))}, and integration by part yields 

Bl ^ ['i/{x){l-'^/{x))dx . (1.41) 


Thus 0 < B^ < oo for any shape ih G 5^. 
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2 Modified Circular Unitary Ensembles 

We present the counterparts of the results of section [L3l for the modified CUEs fdefinition II. 10|1 . Along the 
way we set up definitions and lemmas that will also be used in section |3] and U] Actually Circular ensembles 
can be thought of as simplified models which are helpful to understand the combinatorial structure behind 
the cumulants of linear statistics of the MNS model because no asymptotic estimates are required to pass to 
the limit. In section [231 we review the method introduced in [35]. In section [221 we show that the modified 
CUEs exhibit the same transition as in figure [T1 Finally, in section |2l3l we provide asymptotically equivalent 
kernels for the modified CUEs in the critical regime 5 = a and we deduce a limit theorem from the results 
of section |4| 

2.1 Soshnikov’s method: proof of lemma 11.111 

In [36] lemma I, Soshnikov proved that the cumulants of linear statistics of a determinantal process are given 
by 

|m|=n 

where K is the correlation kernel of the process (the sum is over all compositions m of n). See also equa¬ 
tion (|I.33|) in section Fl.61 for a formal derivation and an explanation of the notations. Applying this formula 
to a modified CUE and using some elementary Fourier analysis we obtain lemma II. Ill 

Proof of lemma U . 1 1[ The correlation kernel of the modified CUE is {x,y) = and for 

feez 

any composition m of n, by formula (11.3411 . 


^(m) _ 

TT[r-K^r^---nK^]= y (2-i) 

KgZnm)i=l 

where by convention kg = For any indices s,r G X and m G N, we know that 

fm(^g_r)= Y f{ki-r)f{k 2 -ki)---f{s-km-i)- 

fcGZ™-i 

Let rfij = mi -I- • • • -I- as in definition (11.351) . For any * = 0, • • • we can write 

J^i+l ^ ^ f + l )f (^mi-t-2 ‘ ‘ ‘ f {krnij^\ —l) 

and if we make the change of variables Ki = k^. in equation (ED), putting everything together we get 

n ^(m) 

Tr[/™^ ■ ■ ■ nK] = E n ) n pL, ■ 

/cGZ'‘+i:fco=fcn i=l i=l 

We can also make the change of variable Ui = ki — ki-i in the previous sum which maps {k G IP : fcg = fc„} 
into {{ko,u) S Z X Zq} and we obtain 

^(m) 

Tr[r^K^r^---nK^] = y E U pZ+u,+-+u^, • 

i ko^'L 2=1 

Hence lemma FT. 11 1 follows directly from formula (11.331) . □ 
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Observe that we recover lemma 1 in [35] for Dyson’s CUE by taking = l|fe|<Ar, 


since then 


i-i 


En Pk-\-ui-\ -j-'u, 




i —0 


2N + 1 — maxjwi + • • • + ujn} 
o<i<r 


+ 


max| — lii — • • . — Urn-} 
o<i<r 


( 2 . 2 ) 


In section [Til we used lemma 11.111 to show that a modified CUE has non-Gaussian fluctuations at the 
macroscopic scale. In the sequel we will use it to investigate fluctuations at mesoscopic scales. Let / G C'o(M), 
0 < 5 < 1, and recall that T = [— 5 , 5 ] with the endpoints identified. When the parameter N is sufficiently 
large, the function f{-N^) is supported in [—\, 5 ] and it can be extended to some function fs £ C'(T). Thus, 
the Fourier coefficients of fs are given by, for any u G Z, 

Mu) = N-^f{uNM . 


Hence, 




H- \-Ur^ 


i 


|m|=n 


k^Z i—1 


(2.3) 


2.2 Central Limit Theorems 

From now on, we will assume that the spectrum of the modified CUE correlation kernel is given by = 
I see (|1.23l) . Moreover, to keep the notations simple, we will write pk instead of p^. 

We start by proving a classical CLT at the Poisson scales {S < a). The proof relies on a simple variance 
computation. Observe that the asymptotic variance of theorem 12.11 matches that of theorem 1 1.31 only up to a 
multiplicative constant. The difference is due to our normalization. Namely, the scaling (11.231) implies that 
the modified CUEs have density 2N at the origin and 

[Efs] = \2N+ O (7V“)| N-^ [ f{x)dx . (2.4) 

P t Af —>00 ) J 

Theorem 2.1. Consider a modified- CUE with correlation kernel (|1.22l - ri.23l) and let f G Hq (R). For 
any scale 0 < <5 < a < 1, the centered and rescaled linear statistic N^~'E.fs converges in distribution to a 
Gaussian random variable with variance 2rH|,||/|||2(R)- 

Proof. When n = 2, equation (IT31) reads 

C],^4Efs] = 7V-2^^/(uiV-^)/(-^A'-')^Pfe(l -Pfe+u) . (2.5) 

u^Z k^Z 

We let al = pk{l — Pk) for any fc > 0. Recall that p-k = Pki then for any n G Z, 

- Pk+u) =Po(l-Pu) + y^Pfc(2 -Pfc+« -Pk-u) 
k^Z k >0 

= po{l -Pu) + 2'^al + '^Pk{2pk -Pk+u -Pk-u) 
k >0 k >0 


Since the coefficients pk G [0,1] and the shape dt is non-increasing, we can check that 


and we get 


^Pk{l-Pk+u) -2j2^k 

k^Z k>0 


< 2\u\ + 1 


'^Pk{Pk -Pk+u) 
0<k 


< |u| 


( 2 . 6 ) 
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If we combine this estimate with formula (12.51) . since the test function / is real-valued, 


<7V-^^|/(nIV-^) (2|n| + l). 


u^'E 


N^Cl. [Efs] - 2N-^ Y. |/(«^■')| E ' 
u^E fc>0 

If we assume that the test function / G the r.h.s. satisfies 

I 2 

N-^ Y \fiuN-^) (2|u| + 1) < CN^ / \f{v)\^\v\dv 

uGZ 

Then dMl) yields 

iV^-“C^.[S/,] = 2IV-^^|/(z,iV-^)' iV-“^a2+ O (iV^-“||/||W) • 


(2.7) 


2V . 


u^E 


k>0 


Moreover, using (11.411) . a Riemann sum approximation gives 

hm N-^Y^I= [ 'P(tT-')(l-4/(tr-i))dt = rB| , (2.8) 

and when S < a, 

/ OO 

-OO 

Since the variance of the random variable Efs is diverging like N°'~^ and its expected value is of order 
by equation (j2.4|l . the CLT follows from Soshnikov’s theorem 1 in [36]. □ 


Note that using the upper-bound (12.6p and the limit (12.81) . we get 


YPk{^ -Pk+u) < 2\u\ 
fcGZ 


O (iV“). 

N—^oo 


Hence, by formula (|2.5|) . 

C^„[S/5]<27V-^^|/(uiV-^)f ||u7V-V O (7V“-^)| . (2.9) 

P \ \ K N^oo ) 

uGE 

/2 

This implies that for any / S Hq , the variance of the linear statistic ^fs remains bounded in the regime 
^ > a. Actually, if (11.231) holds, we have in the regime S > a, 

lim YaTiiNlEfs] = ||/||^i/2 ■ (2.10) 

iV ^oo ^ 

This suggests that at any scale 6 > a we should have the same limit theorem for the modified CUEs as for the 
mesoscopic CUE and the sine process. We can prove formula (I2.10|) in the same way we got proposition 12.II 
but the argument is quite technical and it becomes really sophisticated if we are interested in computing 
the limits of the higher-order cumulants. A better approach consists in deducing the CLT from Soshnikov’s 
theorem [3^ by proving that the cumulants of a given linear statistics have the same limit regardless of the 
shape 'k of the modified CUE. 

Theorem 2.2. Consider a modified CUE with correlation kernel (11.221 — IT.231) and let f G Hq (K). For any 
scale 1 > J > a > 0, the linear statistics 'E.fs converges in distribution to a Gaussian random variable with 
variance ||/||^i/2. 
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Proof. Let us decompose 


( 2 . 11 ) 


Pk — l|fc|<Ar + • 

By assumption, Sk = — 1 when |fc| < N and Sk = when |fc| > N. We can write 

[Efs] = Cl. [Efs] + £lif, 6, a, M-) , 

where £l{f, 5, a, 5') collects all the term which contains at least one factor Sk when we insert the decompo- 

i-i 

sition (12.1111 into equation (12.3|) and expand the products _. Plainly, all other terms exactly 

add up to Cl.[Efs]. Since \pk\, |efc| < 1, we get 


i=0 


\£^{f,S,a,^)\<Cr.N-^^ 








fcGZ 


where Cn = '^j ^ 




. Moreover, by the definition of e^, we have the two estimates 

rO 

\£f^\= Y 1 - < CiV" / 

0<k<N -N<k<0 

poo 

^ |e^,| / «'(t) 

-\T ^ u r\^ u 


dt . 


N<k 


0<k 


Both integrals are finite since ^ G S and there exists a positive constant > 0 such that 


£^U,6,a,'i>)\<C'^N' 


</ ATa — S 






Therefore, all the cumulants Cl. [.ni/a] and Cl. [^fs] have the same limit and the CLT follows from Theo¬ 
rem 1 in ESI. □ 


Remark 2.3. In the terminology of definition 11.81 we have proved that the rescaled correlation kernels 
N~^Kp {N~^x, N~^y) and N~^K^{N~^x,N~^y) are asymptotically equivalent when the condition 6 > a is 
satisfied. We could also have deduced this fact from lemma [2761 below by checking that the CUE kernel , 
see (11.241) . satisfies the property L^B at any scale. 


2.3 The critical regime 

It remains to look at what happens at the critical scale S = a. We have already seen that the variance remains 
bounded as N ^ oo. We can compute its limit by applying a Riemann sum approximation to formula 
For any 0 < a < 1 and any r > 0, 


lim cl. [^fa] = 2 

N^oci ^ 


[ f{u) 

Jr 




1 - T 


t + u 


dtdu . 


( 2 . 12 ) 


Because of some subtle cancelations, it is difficult to use equation (12.3p to compute the limits of the higher- 
order cumulants by Riemann sum approximations. Another approach is to rewrite the correlation kernel of 
the modified CUE before computing the cumulants. From definition 11.101 a summation by parts gives 


Kp {x, y)=Y 

k—0 

1 


{Pk -Pk+i)K’'{x,y) 

k + ^k\ sin {{2N + 2k+ l)Tr{x - y)) 


rW 


OO 

E 

k=-N 




TfV° 


sin(7r(a; — y)) 


(2.13) 
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where G (0,1) by the mean-value theorem. We can use formula (12.1311 to relate the kernel to the sine 
kernel and we will be able to use the ideas of [3S] to compute the limits of the cumulants of linear statistics 
of the modified CUEs. 

Proposition 2.4. At the critical scale 6 = a, the modified CUE kernel and the kernel ^ given by 
(|1.21l) with ri{k) = (TV -|- fc -|- are asymptotically equivalent in the sense of definition \1.8[ 

By proposition 11.91 a similar approximation holds for the modified GUEs. There is only a minor difference 
in the definition of the function rj and the limits of the cumulants of both models will be computed in a 
common framework in section |4l In order to prove proposition 11.91 and 12.41 we need to give a criterion to 
check whether two kernels are asymptotically equivalent. First, we need a new definition. A similar concept 
was introduced in |30) to control cumulants of some complex determinantal processes. 

Definition 2.5. A family of kernels {L^)n>o satisfies the property L^B if there exists a sequence of functions 
Tat : X —>■ R'*' and v > 0 such that for any compact set A C X and all {x, y) € 

\L^{x,y)\ < rAr(a; - y) , 

and 

l|riv||LqA) = ^o^(|iogivr) . 

Lemma 2.6. Two families of kernels {L^)n>o o^nd cife asymptotically equivalent if the family 

{L^)n>o has the property L^B and there exists k > 0 such that for any compact set A C X, 

snp{\L^ix,y)-K^{x,y)\:{x,y)GA^}= p (iV^) . 

^ ^ Af-^oo 

Proof. If we replace , the general term of the cumulant expansion (11.331) is of the form 

Tr = Tr ^ Tr[jV"'' • • • • 

Note that all terms of the last sum contains at least one operator E^. Moreover, we can assume that the 
test function / is supported in a compact set A and that there exists positive constants C and n such that 
sup {I(a;,?/)I : {x,y) G A^} < CN~^. 

Then, if we first consider a trace which contains a single operator , by (11.341) . we get the estimate 

\TTpn ■ ■ ■ J'r‘]\ < CN-^ n ll/ll- / z^)| dp{x^) ■ • • dy{xe) . 

k=2 

Since has the property L^B, there exists T^v : X —)> R+ such that \L^{xk, Xk+i)\ < T^ixk — Xk+i) and a 
change of variables yields 


i 

|Tr[jVi... j7,]| . 

fc =2 

A similar argument shows that each trace which contains j operators E^ is bounded by N~^'^ times a loga¬ 
rithmic correction coming from ||rAr|lLi(A)- Therefore, if we use the notation O introduced at the beginning 
of section \m we get 


Tr ••• = Tr -k O 


(2.14) 


Since the cumulants are linear combinations of such traces, the proof is completed. 


□ 
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Proof of proposition A Taylor expansion of sin(7r(a; 
that 


N-°‘K^{xN-°‘,yN-°‘) 


1 

tN°‘ 


E 


k=-N 


( k + ^k \ 

V tN°‘ ) 


— y)N “) in denominator of formula (I2.13|) shows 


sin {2TTr]{k){x - y)) 
Tr{x - y) 


+ O , (2.15) 

N^oo ^ ' 


and the error term is uniform over any compact subset of Then, by lemma [521 it is enough to prove that 
the r.h.s. family of kernels denoted ^{x,y) satishes the property L^B. Note that the family of kernels 
^{x,y) is translation-invariant on R, so we can just take Tjv = It is well-known that there exists 

a universal constant C > 0 such that for any s > 0 and n > 0, 


r 

sinn(a; — y) 

Is 

(x-y) 


dy < C'log(sn) . 


This implies that 


Since r]{k) = “ -|- 


-N 


< 


c 


- tN' 


E 

k=-N 


k + ^k 

tN°‘ 


log(s77(fc)) 


(2.16) 


pOO 

and / $(t) logt dt < 00 , we deduce from the estimate (j2.16l) that there is a 

Jo 

constant C' which only depends on the shape dt such that ||< C"log(sA^). 


□ 


If we combine proposition 12.41 and corollary 14.101 we get the following limit theorem for linear statistics of 
the modified CUEs at the critical scale. In particular, up to a scaling, it is the same limit theorem as for the 
critical modified GUEs, theorem ll.61 

Theorem 2.7. Let f € id^/^(M) with compact support, 0 < a < 1 and 4/ G jj'. The linear statistic Sfa 
of the determinantal process with correlation kernel (11.221 — Tl. 231) converges in distribution as N ^ oo to a 
random variable ^^.rf whose cumulants are given by 


C" [E^.rf] =2 tBI [ fitrdt 

Jr 



n./(M*)$(a:*) 


G’^{u,x)d^-\d^x . 


3 Central Limit Theorems for the Modified GUEs 

We begin in section 13.11 by proving some technical lemmas that are needed to compute the asymptotic 
variance of linear statistics of the modified GUEs. In particular we get formula (11.151) for the variance at 
Poisson scales. In section l5151 we prove theorem ll.4l bv comparing the rescaled correlation kernel of a modified 
GUE to the GUE kernel using the perturbative method developed in section (2731 All these results are based 
on the asymptotics of the Hermite functiions and the GUE kernel which are presented in section |5l 

3.1 Proof of theorem 11.31 

We start by proving a classical formula for the variance of linear statistics which is valid in a general context. 

Lemma 3.1. Given a determinantal process with a correlation kernel K of type (HU, for any test function 
f G C'o(X), we have 


Var^ [S/] 


E'^fe / ff {f{x)-f{y)f\K{x,y)\'^n{dx)n{dy) , 


where al = pk{l - Pk)- 


19 

























Proof. If we apply formula (11.331) when n = 2, 


Varx [S/]=y f{x)‘^K{x,x)n{dx) - JJ f{x)f{y)K{x,y)K{y,x)n{dx)n{dy) 

= \ JJ (fix) - f{y)fK{x,y)K{y,x)y{dx)n{dy) 

- If f{x)‘^K{x,y)K{y,x)y{dx)y{dy)+ J f {x)^K{x, x)y{dx) . 

Note that when the kernel K is reproducing, the last two terms cancel. In general, since the function ifk are 
orthonormal, we get 


JJ f{xfK{x,y)K{y,x)y,{dx)y,{dy)+ J f {xf K{x, x)y{dx) 

■--'^PkPj J (pj{y)T^{y)p{dy) J f{xfipk{x)^{x)y{dx)+'^pk J f{xfipk{x)i^[x)y{dx) 

k,j k 

■^<^1 f fixf\pkix)\^Pidx) . 


□ 


For linear statistics of the modified GUEs, there is no counterpart of equation (12.31) . but we can use lemma [01 
to compute the asymptotic of the variance. We call the reproducing variance the quantity 

Mf) = \ JJ \f(.x) - f{y)\^\K^^^\^ dxdy . (3.1) 


This definition comes from the fact that, if the correlation kernel K is reproducing, then = 0 for all 
k G Z_|_ and Vari^ [S/] = Vb(/) for any linear statistic. On the other hand, we call the Poisson variance 
the quantity 

OO 

VAf) = T.'^k 

k=0 

This definition is motivated by the observation that considering a constant test function, we get Vari^ [ff\ = 
Vo-(l) and this quantity measures the extra randomness induced in the process from the fact that the corre¬ 
lation kernel is non-reproducing. In particular for a modified GUE we have 


f{x)'^\(pk{x)\'^dx 


(3.2) 




k-N 

tN°‘ 


1 - ^ 


k-N 

tN°‘ 


(3.3) 


and consequently by dSS, 


Var^^^[#] - TiV“ / «'(t)(l - ^{t))dt = TiV“ B| . 

J —OO 

This is equation (11.101) . see also equation (11.91) in the introduction for a probabilistic interpretation. We shall 
see that, except at the critical scale 5 = a, only one component of the variance is asymptotically relevant. 
We begin by computing an asymptotic formula for the Poisson variance. 

Lemma 3.2. For any {a,S) G (0,1)^ and for any function f G Co(U) we have 

VMs) = [ f{xfdx+ o (1V“-^) . 

2 AT-i-oo 
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Remark 3.3. It is not difficult to adapt the proof so that Zemmo[321 holds for any function f € L^(]R) which 
is uniformly continuous. In particular, by Morrey’s inequality, this covers all test functions in the Sobolev 
space 


Proof. Recall from definition 11.11 that ipk{x) 



are the rescaled Hermite functions and 


ll</5fc|lL2 = 1. Since we assume that the test function / is bounded, by equation (13.21) . for any 0 < e < 1, 


Vaifs)= <^N+k f fs{xf\<PN+kix)\'^d.x + 0 


I oo ^ 

|fc|>Ari- 


<XN+k 


The condition vb G 5^ guarantees that the error term is converging to 0 as iV —>■ oo. Actually, under the 
stronger assumption 'b G 5^', this term decays faster than any power of N and it will be neglected in the 
sequel. Moreover the assumption that / has compact support in conjunction with the condition |A:| < 
implies that we can use the bulk asymptotic for the Hermite functions, equation (15.51) . 

Namely for any x G supp(/). 


ipN+k{xN ^) = cos\{N + k) - F{xnN '^))j + O {N ^) , 

L \ Z / J n—>-oo ' ^ 


(3.4) 


where we set Then 


VMs)=N-^ ^ 


^N+k 


f{xf 


|fc|<Ari- 


cos 


(iV + fc)(|-F(a;7vA^-^))]pdx + 0 ^ 


'N+k 


|fc|<Ari 


Observe that according to (13.31) and (11.101) . ^ ^‘N+k — Var^w [ff\ = 0{N°‘), and the previous estimates 

gives 


vAfs) = 


N 


-s 


Var^« J#] ll/ll i.+ I fixrcos[2{N + k)F{xNN-A]dx'^+0{N^-^A- 


|fe|<Ari- 


The second term is a sum of oscillatory integrals and we will show that it converge to 0 as —>■ oo. Let us 

make the change of variable z = N^F{xnN~^)- By definition 15. 11 


f (x)^ cos [2{N + k)F{xNN ^)] dx = — 




^1 + G' (zN-^) cos [2(7V + k)N-A] dz 


Since the function / is uniformly continuous (we assume that / has compact support in there exists 

a sequence cat \ 0 such that uniformly over all |z| < L and all |fc| < 

'2N^ 




< cn ■ 


Since G"(0) = 1/2, it follows that for any |fc| < 

J f{xA cos [2(A^ + k)F{xNN~^)] dx = — J ffi ^ cos [2(A + k)N~^z^ dz + O {e^) . 

Since the sequence {N — —>■ oo as —>■ oo, by the Riemann-Lebesgue lemma, we can also assume 

that 

sup [ f{xA cos [2{N + k)F{xNN~^)] dx < cn . 
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Going back to the Poisson variance, we have shown that 

VAfs) = ^Var^^ J#] {||/||2, +^o^(e^)} 

The lemma follows after replacing Var^w [#] using equation (11.101) . □ 

In order to prove (11.151) it remains to estimate the reproducing variance Vo(/). Before proceeding we need to 
recall a few properties of the GUE correlation kernel dH. We refer to section [5] for more details, although 
the normalization are different. We deduce from the Hermite functions uniform bound (15.41) that there exists 
a positive constant C such that 

II^^mIIoo < . (3.5) 

This yields a trivial uniform bound for the GUE kernel 

M-l 

ll^(f L ^ E < CMN^/^ . (3.6) 

The connection with the modified GUE kernel comes form a summation by parts, 


= Y^(Pk - Pk+i)KAx,y) 

fc=0 


1 

tN°‘ 


E ^ 


k=-N 



r^N+k 


ix,y) . 


where $ = —'k' and G (0,1) are given by the mean-value theorem. If we further suppose that the shape 
'k S using the uniform bound (13.61) . we can show that for any P > 0, 

K^^Ax.y)=r-^N-‘^ Y, + . (3.7) 

|fc|<riV“ ^ ' °° 

Remark 3.4. The assumption d)* S g' implies that, choosing T = (logA'^)^, the error in formula (13.71) decays 
faster than any power of N. This approximations hold for more general shapes, though not necessarily at all 
mesoscopic scales with such a good error term. Moreover, the condition 'k S g' makes the proof quite simple, 
otherwise we would need to take into account the speed of decay of and make more precise estimates. 

Lemma 3.5. For any 0 < a < 1 and any scale 0 < d < 1, there exists a constant C > 0 such that for any 
1 /2 

function f S Hq' (R), the reproducing variance satisfies for any N > 0, 

Voifs) < C\\f\\Y, . 

Proof. To simplify the notation, let us assume that the temperature r = 1. We can also neglect the parameters 
fk = 0 and the error term in formula (13.71) . Indeed, taking e.g. P = (loglV)^, this error decays faster than 
any power of N. We will denote xn = xN~^, resp. yN = yN~^. Let us fix a parameter L so that 
supp(/) C [—■§, ■§]• According to the sine-kernel approximation (15.101) . if the density N is sufficiently large 
compared to L, there exists a constant C > 0 such that 

\N-^K^+^{xN,yN)\<C\x-y\-^ 

for all x,y G [—L,L] and all |A:| < rA^“. Since 

If 

|fc|<rAr“ ^ 
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this implies that 


< iV-“ ^ $ (fc7V-“) \N-^K^+>^{xN,yN)\ < 

|fc|<rjv“ 

Hence, by formula (EU and a change of variables, 


C 

\x-y\ 


yo{fs) = \ JJ \fs{x)-fsiy)f\K^,aix,y)\^dxdy 





fix) - f{y) 

x-y 


dxdy+ JJ \f{x) - f{y)\'^ \N ^K^^^{xN,yN)\^ dxdy . 

R2\[-L,L]2 


(3.9) 


For any L > 0 the first integral in (13.91) is bounded by ||/||^i /2 and it remains to show that the second 
integral can be made arbitrary small. We claim that there exists a constant C{f) > 0 which only depends 
on the test function / such that 


JJ \fix)-f{y)\^\N-^KlJxN,yN)\"dxdy<Cif)L-^ . 

R2\[-L,i]2 


(3.10) 


Thus, choosing the parameter L sufficiently large in (13.9p . this implies that the variance Vb(/ 5 ) is bounded 
by Cll/ll^i/a. The rest of the proof is rather technical, it is devoted to the estimate (j3.10|) . First, we see that 
by equation (lT7l) . 


Ko~^''ix,y)KQ^^{x,y) 

\k\<rN°‘ 

bl<rAf“ 


i^^,«(x,y)|'<fV-2“ ^ $(fc7V-“)<i>(jiV-“) 


Hence, since the function / is supported in ■§■] and by symmetry of the kernel we obtain 


JJ \fix) - fiy)f \N ^K^^^{xN,yN)\^ dxdy 


(3.11) 


-L,L]2 


< $(fcAf-“)$(jW-“) 


|fc|<rAr° 

bl<rw“ 


fix) - fiy) 


\v\<L/2 

L<.x 


X-y 


iXN — yN)‘^Ko~^'‘ixN,yN)Ko~''\xN,yN) 


N+j, 


dxdy . 


The bulk asymptotic (15.21) implies that there exists a universal constant C > 0 such that for all |fc| < FiV“ 
and for all |a;| < 

\^N+kix)\ < J^2)1/4 (3.12) 

In particular, for any \y\ < L/2, we have \ipN+kiyN)\ < C, and by the Christoffel-Darboux formula (15.61) . for 
any a; G R, 


|(xAr — yN)'^Ko~^’^ixN,yN)Ko~^^ixN,yN)\ 

< (\ipN+kixN)‘PN+jiXN)\ + \<PN+kiXN)‘PN+j-lixN)\ + \‘PN+k-lixN)y’N+jixN)\ 


(3.13) 

IV^Af+fe-l (xAr)v3Af+i-l (xAr)l) 


Dehne 


Jk,jiL) 



\y\<L/2 

L<.x 


fix) - fjy) 

x-y 


2 

\ipN+kixN)‘PN+jixN)\ dxdy . 
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According the estimates (13.8L (I3.11|) and (13. 1311 . we see that in order to prove (j3.10ll . it is enough to show that 
there exist a constant C{f) such that, for any |fc|, \j\ < rN°‘, the integral Jkj{L) are bounded by C{f)L~^. 
Since supp(/) C [—-j, -j], we have for any |a;| > L, 


fix) - f{y) 

x-y 


< 1 


yGsupp(/) 


2||/l|oo 

\x-L/2\ 


< 1 


i/Gsupp(/) 



and, if C'(/) = 4||/|1^| supp/|, we get 


rL/2 

fix) - fiy) 

J-L/2 

X-y 


dy < 


Cjf) 

\rp\2 


Hence, by a change of variables, 


JkAL) < Cif)N 


-s 


! LN-^ 


\ipN+k{x)<fN+jix)\ 


dx 


(3.14) 


(3.15) 


We can use the asymptotic formulae of section [S] to estimate this integral. At this point, we need asymptotic 
expansions which are valid up to the edge. In the bulk, we can use the upper-bound (13.121) . On the other 
hand, by formula (15.31) and since < H{x), we deduce that there is a constant 0 < /3 < such that 

for all |fc| < TN°‘ and for all x > 


\‘PN+kix)\ < . 


Hence, it follows from (13.151) that 
JkAL) <Cif)N-^ ' 


dx 


'LN- 


^\/4 — TT^a;^ 


ii'^iv+fciiooiiv’iv+.iioo / ' r 

I 2 ( 1 -N-^/^) X‘^ I 2(l + iV-l/ 2 ) 


The third integral is bounded by e , By (|3.5I) . ||</?Af-i-fe||oo < and the second term is of order 

7V1/3-1/2 _ TV” 1/6^ xhe first term gives the leading contribution, namely 


2(1-N 


-iV-V2i 


N 


-s 


dx 


< LN- 


2\/4^ 


zdx < 


C 


We conclude that JkjiL) < C{f)L ^ uniformly over all |A:|, |j| < rA^“ and the estimate (13.101) follows by 
combining (|3.8I) . (13.111) and (13.131) . □ 


We can now conclude the proof of formula (11.151) and hence of theorem 11.31 It follows immediately from 
lemmas ITTI 13.21 and 13.51 that in the regime d < a, for any test function f G Hq (R), 

Var^^ = / f{x)Ax+ o (iV“-^) . 

2 Jr Af-i-oo 

The same argument shows that, in the regime 6 > a, 

Var^^ jS/,] < C(||/||i. + ll/ll?,,/.) (3.16) 

At the GUE scales (<5 > a), the limit of the variance is given by theorem 11.41 which is proved in the next 
section. At the critical scale d = a, by lemma 1?^ the Poisson variance 14(/a) converges to ^ B^ 11/11^2 and 
the limit of the reproducing variance Voifa) is computed in appendix lAl bv a Riemann sum approximation. 
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3.2 Proof of theorem 11.41 

Theorem 3.6. Let S be the GUE eigenvalue process with correlation kernel Kq given by m- For any 

1/2 

0 < 5 < 1 and any function f £ Hq' (R), as the number of eigenvalues N ^ oo, 

Efs A/" (0,11/11^1/2) . 

Theorem 13.61 was first established in for the resolvent function x (x — z) ^ where Sz > 0. A general 
proof was given only recently in [20] for test functions with compactly supported Fourier transform. Their 
proof exploits a nice connection between the characteristic polynomial of a GUE matrix and a log-correlated 
Gaussian process. In [4], a generalization of theorem l3.6l is given for Gaussian /3-Ensembles and test functions 
/ S C'o(R). When the sampling scale 0 < <5 < 1/3, the analogue of theorem 13.61 holds for Wigner matrices 
as well, see m- In [S] , yet another analogue of theorem 13.61 was established for a large class of Orthogonal 
Polynomial ensembles and test functions / £ (70(1^). Their results covers the Jacobi Unitary Ensembles 
but not the Gaussian Unitary Ensembles. In the paper [26] in preparation, we give an elementary proof of 
theorem 13.61 based on the cumulant computations presented in section |4l Actually, our method applies for 
any test function / £ C'g(]R) and for any determinantal process whose correlation kernel satisfies the sine- 
kernel asymptotics of theorem 15.21 In particular, based on the results of [25], this includes the Orthogonal 
Polynomial ensembles discussed in |8] . 

We turn to the approximation of the modified GUE correlation kernels at the GUE scales {6 > a). By 
definition 11.81 proposition 13.71 below and theorem 13.61 imply the central limit theorem 11.41 

Proposition 3.7. For any shape dt £ , the modified GUE correlation kernel ^ and the GUE kernel 

Kq are asymptotically equivalent at any scale S > a, i.e. 

N-^K^^^{N-^x,N-^y) ^ N-^K^{N-^x,N-^y) . 

Proof. We can assume that r = 1 and we fix a parameter L > 0. The condition dt G 5^' implies that 

^{kN-°^)+ Y (1 - . 

fc>rA^“ k<-rN^ 

So that, if we choose T = (log A^)^, these sums decay faster than any power of N. This fact combined to the 
uniform bound ()3.5p implies that 

N-l FN^ TN°‘ 

- Y Pk{x)(pk{y) + X! ^ {kN~°') ipN+k{x)(pN+k{y) - X! ^ {-kN~°‘)) ipN-k{x)(pN-kiy) 

k=0 k=0 k=l 

with a uniform error of order O that we will neglect. If the parameter N is sufficiently large, 

it follows from the bulk estimate (|3.12l) that all x,y £ [—L,L], 

Y \pN+k{N-^x)ipN+kiN-%)\ < 2C^rN‘^ . 

|fc|<rAf“ 

Since G [0,1], using the notation O introduced in section [L6l this implies that 

Af-l 

N-^K^^{N-^x,N-^y) = N-^Y‘l^k{N-^x)MN-^y)+ O (iV“-^) 

k=0 

uniformly for all x,y £ [—L,L]. The sum in the r.h.s. is the rescaled GUE correlation kernel (11.411 . By 
lemma 1^61 to prove that the kernels ^ and Kq are asymptotically equivalent, it remains to show the 
latter satisfies the property L^B. Taking M = N in the approximation (15.1011 implies that there exists a 
positive constant Cl such that for any 0 < J < 1, 
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N-^K^{xN-^,yN-^) \ < r^(x- y) = Cl 


'n if lx - y| < 


1 

\x-y\ 


We can immediately check that J TN{z)dz < 4 Cl logiV and the proof is complete. 


□ 


4 Cumulants of the Critical models 

In this section, we prove theorem ll.6l and l2.7l then we analyze the random processes which arise from the 
critical modified ensembles. Actually we will not investigate directly the modified ensembles but the process 
with kernel given by (11.211) . By Dror)ositions ll.9l and l2.4l there are two particular choices of the function 
r] which correspond to the modified GUEs, resp. CUEs, but our analysis works as long as rj satisfies the 
conditions (14.41) - (14.51) . In section im we show that ^ is the correlation kernel of a determinantal process 
and we prove proposition [Tj^l The convergence of smooth linear statistics of these processes is established in 
sectionsee corollary 14. 101 In section H751 it is proved that the random variables 7 ./ are not Gaussian 
if r > 0. Our argument holds for almost all shapes 'k G S'. But, unfortunately, it does not apply to the MNS 
model because of its special property; see proposition 11.71 In section ITTl we show that the MNS ensemble 
at the critical scale, despite this property, converges to a random process which is not Gaussian. 


4.1 Asymptotically equivalent kernels for the critical modified GUEs 

Lemma 4.1. Let N,t,T > 0, a G (0,1), 4' G 1? and y be a non-decreasing function. The kernel L^^ given 
by (ll. 2 ip defines a translation-invariant determinantal process on K. 

Proof. The fundamental property of the kernel L^ ^ is that it is translation-invariant. Hence we can define 
its Fourier transform 

X! ( .j-jya ) ■ (^-1) 

|fc|<rAr“ ^ ' 

Plainly, the function G L^(IR) and by (11.211) . 

L^.r,i^,y)= f . (4.2) 

Jm. 


This definition was used in |36j where it is given that the condition for any translation-invariant kernel L to 
define a point process is that its Fourier transform satisfies 0 < L < 1. The parameters have been chosen 
so that 7 ]^$ () = 'k ( 7 ^) — 'k ( 7 ^) and it follows that for any u G K, 


L^Jv) < 4-(-rT-i) < 1 


Moreover, since $ > 0 by assumption, ^ > 0 and we conclude that L^ ^ is the correlation kernel of some 
determinantal process. □ 


Proof of vrovosition \l.fh Let 1/3 < a < 1, 4^ G S'', and T = (log A^)^. We also assume that r = 1 to simplify 
the notation. We combine the approximation (ITTI) of the modified GUE kernel ^ with the asymptotic 
formula of theorem 15.31 with xq = 0. Namely, taking d = a in formula (15.111) . we get 


N-^K^^{xN- 


,yiV-“) = fV- 


^ E 

|fe|<rAf“ 


$ 


fc + 6 \ sin [tt“ y'l -b k/N[x - y)] 


N° 


7r(x - y) 


-b O (a^i-3“) 
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The estimates (j2.16|) shows that, with r]{k) = + k/N, the kernel given by (11.2111 satisfies the 

property L^B and it follows from lemma [27^ that 

□ 

Proposition 11.91 and 12.41 imply that each of the modified ensembles have an asymptotically equivalent kernel 
at the critical scale of the form ^ with 

i) r]{k) = (N + k)N~°‘ for the modified CUEs , , . 

a) ri{k) = i_/v^““y^l + k/N for the modified GUEs . 

In the sequel, we will compute the limits of the cumulants for any determinantal process with kernel ^ 
which satisfies the following conditions. The function r] is non-decreasing and it satisfies uniformly for all 
|fc| < TN°^, 

y{k) =N'' + pkN-°^ + O , (4.4) 

N^oo 

where z/, /3, e > 0 such that N'^ T and 

lim 7V‘'max{4'(r),l - 4'(-r)} = 0 . (4.5) 

N^oo 

In particular, for the modified GUEs (resp. GUEs), the asymptotic expansion (14.41) holds with v = 1 — a and 
/3 = 1/4 (resp. /? = 1) and, if the shape 'h G i?', the condition (14.51) holds for any a G (0,1) with T = (log A^)^. 


4.2 Proof of theorem 11.61 

Given the expression (14.21) of the correlation kernel we can repeat the proof of lemma fl . 1 II replacing 

sums by integrals and we get the following formula 








^ / II X! M(™) / H - \-Urni)dv\(P 


(4.6) 


|m|=n 




where the sum is over all compositions m of the number n\ see (11.321) . Combining this formula with (14.11) . 
we get another expression for the cumulants that will be appropriate in the limit IV —^ oo. In this section, to 
simplify the notations, we will assume that = 0 and (unless stated otherwise) all sums run over |fci| < riV“. 
Define for any composition m of n > 2, the function 


Hm (w, k^ 



|ll-|-lll-|- 


dv . 


(4.7) 


Lemma 4.2. For any function f G C'o(K.), the cumulants of linear statistics of the determinantal process 
with correlation kernel (11.211) are given by 



where the function Hm is given by (14.71) . 

Proof. To simplify the notations, let us assume that r = 1. 

By formula 611), for any composition m of n of length i and any n G we have 


2=1 


= N 


— Oi^ 


E 

i,-" 2=1 


kj 




(4.9) 
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If we let = 0 and for any 1 < i < n, 


,kn ^ 


we get 


E n*(jd;)n* 


2=1 


ki,‘“ ,fcn ^=1 


ki 




2=1 


By (14.7L this implies that 


f ^ ^ / L \ 

•'**=1 ki,--,k„ i=l ^ 

r ^ 

Observe that by (14.9L since iV““ ^ ^{kN~°') < 1 and condition (14.4p holds, / L^,^{vi)dvi < AH'". 

K. •'* i=l 

Moreover, by definition 0 < < 1 — 4>(—F) and according to condition (14.51) . 


lim 

N^oo 


+ Ml + • • • + UrrTi )dv — 0 . 


Thus, by formula (I4.10F 


i=l 


ki<.‘“<.kn 2=1 


,j(m + Ml + • • • + Mmjdw ~ TV ^ ^ Hm(M, crfc). 

i7GS(n) 


We conclude by using formula (14.61) . 


□ 


We can use the notations (11.361) and (11.371) to compute the function Hm given by (14.7F The computation 
is analogous to the proof of formula (j2.2l) . In the sequel, we will always use the conventions £ = £(m), 
= A™(m) and s =^(„)(cr). 

Lemma 4.3. Let k G Z” and u G Rq. For any a G S(n) and any composition m of n>2, 

+ 


Hm(M, Crfc) = 


‘^■n{ka{s)) - max - r?(fc^(,)) + r]{k^(^s))'^ - max - ? 7 (fc^(,)) + 77(fc,^(s))} 


(4.11) 


Proof. Let Vi = Uj. The change of variable w = v — Vs^ in (I4J1) gives 


i=i 


Hm(M,crfc)=/ TTl|u,+A”|<r)(fe„(i)) rfw; = H u; : |m; + A™| < r7(A:<^(i)) 

•'* i=l i=l ^ 

By definition mini<^{cr(i)} = cr(s), and since the function ry is non-decreasing, for any k G Z", 

V{ka(s)) = min77(fc<,(i)) . 

Then, since A™ = 0 by (jl.36|) . we get 

^ r 

Pi < I MI -I- A™| < ri{kcr(i)) 


2=1 


(4.12) 


(4.13) 


-dik^is)) + max { - A™ - r]{k^(^i)) + vik^(s))}, dikais)) - max {A™ - 77(A:„p)) -b d{kais))} 

l<l ~ l<l ~ 
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This interval is non-empty if the condition 

> max {A™ - 77(fc^(i)) -f v{ka{s))} + max { - A™ - r]{k^(^i)) + v{ka{s))} 

is satisfied, in which case equation (14.111) follows from equation (14.121) . □ 


We are now ready to prove our main result, i.e. to compute the limits of the cumulants of the linear statistics 
of the modified ensembles by applying a Riemann sum approximation to equation (li^ . The argument is 
similar to the proof of Lemma 2 in |35j but it is more complicated. To keep the proof as transparent as 
possible, it relies on three lemmas which will be proved afterwards. 

Theorem 4.4. Assume that the conditions (lOl) and gSl) are satisfied and let f G Hq{W). For any n>2, 


lim 

N—¥oo 


C” 


N 
^ ,77 


[S/] = 2/3rB^ 


f{trdt -2 



M(m) G™ (u, cc) d” ^udFx 


where the function G^{u,x) and constant are defined by equations (11.381) and (11.391) . 

Proof. Throughout the proof, we will use the familiar inequality (13.8p without any reference. Let |u|i = 
|mi| + ■ ■ • + l^nl and 




/ 1 


V tN‘^ 


E 

fci<...<fc„ i=i 



If the parameter N is sufficiently large, we claim that for any 
that |/c|oo < A^“r, 


(^Vikcris)) - max jA*^ - + ?7(fc,T(s))}^ 

(4.14) 

M G Rg, any a G S(n), and for all fc G Z" such 


H. 


< 


fiu,ak) - 2??(fc„(^)) -brnaxlA)^ - + v{ka{s))'\ -binax|-A'^ - + ?7(M£))} 


0 


if |u|i < ^ 


18|m|i else 


(4.15) 


First note that, since A™ = 0, 



Moreover, by (14.131) . 

- Vika{i)) + Vikais))^ -?<f • 

By formula (I1.36L for any composition m of n, we have |A™| < |u|i for all i < £. Hence, we conclude that 


0 < max I ± A™ 
“ i<e I - 


v{ka{i)) + v{ka{s))'\ 


< lwll 


(4.16) 


When the parameter N is large, condition (14.41) implies that for any |/c| < TA^", 

— < v{n) < 2N'' . 

Thus, if we also suppose that |n|i < by (14.161) . 


Vika{s)) > max jiA*^ - riik^^i)) + vika(s))] 


(4.17) 
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By lemma 1131 we conclude that when litli < 

Hm(u,crfc) - 27?(fc^(^)) + m|x|A'^ - ^7(^0) + ’•'“If “ vika{i)) + viK{s))'\ = 0 . 

For the second estimate, we observe that by formula (14.1111 and the estimates (I4.16|) and (I4.17L 


0 < Hm(w, crk) < 2r]{k„(s)) < 4A^'^ • 

Then, by the triangle inequality and (14.161) . the l.h.s. of (14.151) is bounded by %N'^ + 2|u|i. Thus, we have also 
proved (14.151) in the case when |u|i > If we combine this estimate with formula (14.81) for the cumulants 
of the random variable S/, there exists a positive constant which only depends on n such that if the 
parameter N is sufficiently large. 


+ <Cr, [ (4.18) 

. Jjjn U ' 2 J ^1 

where the function is given by (14.141) . Taking /^ =••• = /„ = / in lemma lT6l below implies that the 

r.h.s. of (j4.18|l converges to 0 as —>■ 00 . Thus the limits of the cumulants is given by 


hm [S/]= hm [ l[f{u.){Tliu) + Tl{-u)}d-^ 

1—^00 ,11 N—¥00 J'^n 

= 2 lim f 5R J n f{u,) \ Tl{u) d^-^u . 


(4.19) 


The next step is to compute the limit of T'^{u) as fV —>■ 00 ; see equation (14.2211 . Observe that, according to 
condition (j4.4ll and since the max function is Lipschitz continuous, we get uniformly for all u £ Rq , 


^lin) = 




E 
E 

(TGS(n) 


h 

^ kcrjs) - kl 


M(m)x (4.20) 

n 

- max I^ ~ 1) + O 

i<l ( ^ ATa j y jv-s-oo '' ’ 


By lemma H75] below. ^ M(m) = 0 and, since they are independent of m, we can remove the two terms rj{ki) 
and kiN~°‘ from formula (14.201) . Hence we have proved that 


'^Niu) = 


( 1 


\tN° 


k\<---<kn 1 ^ ^ |m|=n y (TGS(n) J 

(4.21) 


where G|^-q is given by (11.381) . Then, a Riemann sum approximation implies that 


:= Tf'So(w) = / f[^iz,T ^) M(m)j/3 ^ z^(^) - G™(m, z) j d 

JR'^ Z_1 I_I_ \ / 


|m|=n 


(TGS(n) 


The first term of this equation is independent of the Fourier variable u £ Rq and it can be computed explicitly; 
see lemma 14771 below. Furthermore, making the change of variables Xi = ZiT, we obtain 


„ n 

T^{u)=PtBI- / E M(m)G^,(w,x)d 

i=l |m|=n 


X . 


(4.22) 
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Now, we can deduce the limits of the cumulants of the random variable from equations (14.1911 and (I4.22I1 . 
By (11.381) and the estimate |A™| < |m|i, we get 

sup { G™(a:, u) : x € r > O} < n!|M|i . (4.23) 

Moreover, since |fccr(s)| < |fci| + \kn\ for any k G Z”, by formula (I4.21|l . there exists a constant C which only 
depends on n such that for any > 0, 


P n 

T%{u)\ <c TT $(a:iT“^)(l + \xi\ + \Xn\ + |M|l)d" 


i=l 


(4.24) 


The assumption G guarantees that the r.h.s. of (14.241) is finite. From lemma |T6] below, formula (14.221) . 
and the dominated convergence theorem, we conclude that 


lim Clr, [Ef] = 2 f 5R I n /(u,) I f ff $(x,) V M(m) G^iu, x) ] d^-^d^x 

JK IV j V J^<t\ Ir^n J 

The final observation is that the integral over Rg can be written as a convolution. Namely a change of 
variables gives 


n « n—1 

T[f{Ut)d'^~^U= / f{vi)T\ f{Vi- Vi-i)f{-Vn-l)d'^~^V = f * 


/(O) 


If we replace /*•••*/ = 
theorem 14.41 is completed. 


/" and evaluate at 0, we get 


/R" 


n fiui)d: 


n—1^ 


f{t)^dt and the proof of 

□ 


Now we prove the lemmas that we used to get theorem 14.41 The first lemma is classical, it was already used 
in |35j and also in the context of invariant ensembles in the complex plane [1]. 

Lemma 4.5. For any n > 1, 

M(m) = i5i(n) , and |M(m)| < n!2"“^ . 

|m|=n |m|=n 


We have assumed that our test function / has compact support since the original problem is to study 
mesoscopic linear statistics. However this assumption is not necessary to prove theorem 14.41 We shall 
certainly require that / G and, according to the estimate (I4.18L the regularity condition needed to 

prove theorem 14.41 is that, for any n>2, the integral 


[ f{Un) (l+|M|l)d” < OO 

Jri^ 


A sufficient condition is provided by the following lemma (since by assumption / G L^(R) and ||/||oo < oo). 
Lemma 4.6. Let n > 2, for any functions /i, ..., /„ G 


P n 

/ fi{ui)---fn{un) (1 + |mi|H- \un\)d'^~^u < n 2 ^ Y[ {Wfj\\o° + WfjW' ( 4 - 25 ) 

•'*0 j=l ^ 
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Proof. By the Cauchy-Schwartz inequality, 


[ fl{ui)f 2 ( -^Uj\ (1 + |Ml|)rfui < 2 [ 2II/1II00II/2II00 + /" fl{u)f 2 ( - 


u|fiu 


< 2 


(2||/l||oo||/2||oo + ||/l|ki||/2||L^ 


A similar argument shows that for any 1 < j < n, 

ll/^llLiVlI/.lli. <2(||/,|U + ||/,||^x) . 

Hence, it follows from (14.261) that 


(4.26) 


(4.27) 


[ fl{ui)f2( 
Jr V 


A + |ui|)dui < 4(||/i||oo + ||/l||_f/i)(||/2||oo + ||/2||ffi) 


and, if we combine this estimate with (14.271) . 

[ fl{ui) ■ ■ ■ fniUn) (1 + |Ml|)d””^U < 4(||/i||oo + ||/l || ) (||/2 || oo + IIMIffi) 

yjiloo + ll/jllffi') • 


Li 


j>2 


< 2”]^ 

3=1 


The upper-bound (I4.25|) follows by symmetry. 


□ 


The next lemma shows how the shape-dependent constant defined by (11.391) arises in equation (14.221) . 
Lemma 4.7. For any n > 1, 


B n _ 

— 


(TGS(n) |m|=n 


(4.28) 




n—1 p 

Y^bl / z$(z)vI/(z)'=(l-4/(z))"-'-'=dz , 

k =0 


where Xcr(s) = min {XaU)} by p.37p for any x € K", and the coefficients are given by equation (11.401) . 

Proof. Let P„ be the uniform measure over the group S(n), so that we can view (t(s;) = mini<;{cr(i)} as a 
random variable. Then, we can rewrite equation (14.281) as 


P n 

B^ = n! / V M(m)E„ 

7r" 


X . 


i=l |m|=n 

We claim that for any ^ = 1,...,n and for any k = 0,... ,n — 1, 

''k + l'\l!(n-l)! 


min(T(z) > n — k 


i<i 


I 


n\ 


(4.29) 


(4.30) 


To see this, observe that if A: -I-1 < Z, since there are only k elements in {1,..., n} greater than n — k, one of 
the I first elements of a has to be less than n — k and therefore the probability in question is 0. 

On the other hand, if Z < fc -|- 1, then n — k is smaller than the minimum of the Z hrst entries of cr if and only 
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if these entries are drawn from the set {n — fc,..., n}. Since the order of these entries and that of the {n — 1) 

/fg _l_ i\ 

last entries is irrelevant, the number of such permutations is ( ^ j/!(n — 1)1. 

Hence, by definition and equation (I4.30L the distribution of cr(s;) is given by 

-1 


[cr(s;) = n-k] = 


k 

I - 1 


Then, by definition of M, (jl.32|) . 


M(m)E„ 

|m|=n 


(- 1 ) 


l-\-l 


1^1 


" i: 

|m|=n ^ ^ ^ 

^{m}—^ 


(4.31) 


If we integrate successively over Xi, • • • , Xn-k-i and over x„, • • • , Xn-k+h and use the relationship $ = —4'', 
we find that for any fc = 0,..., n — 1, 


]G[ 4>(xi)x„_fcd"x = 


1 


k\{n — k — 1) 


- [ 4'(x)'=(l-4'(x))"-'=-i$(x)xdx . (4.32) 

■ Jr 


Then, if we combine equations fl4.29|) . (|4.3ip and fl4.32p . we get 






ni n 




l\ \ l 


-1 


n—1 


i: JEwv 


{k — l + l)!(n — k — 1) 


T [ 4'(x)'=(l-4'(x))"-'=-i4>(x)xdx . 

■ Jr 


n)=/ 

, -1 


We see that we can simplify ^ (") from the previous formula and exchange the sums over I and k. In the 
end, we obtain 


n—1 


=i: 


A;=0 ^ Z=1 


\Z+1 


n — I 
k + 1 — I 




I m| —n 


If we define the array according to (jl.40l) . the lemma is proved. 


□ 


The Cumulant problem is generally not discussed directly in the literature, so we provide a simple criterion 
which guarantees uniqueness of the law of a random variable given its cumulants. 

Lemma 4.8. Given a sequence of random variables X]\j whose Laplace transform is well-defined such that 
for any n > 1, the cumulant —>■ as N ^ oo. If there exists constants c,v > 0 such that 

I I < cn!x" , 

then there exists a random variable X^o whose cumulants satisfy C"'[Xoo] = and the sequence X^ ^ X. 


The condition of lemma 14.81 is very natural and its proof follows from a straightforward adaptation of the 
argument that is used when dealing with the Hamburger moment problem (see e.g. section 3.3.3 in |15jl. 
Next, we use this criterion to deduce from theorem 14.41 the weak convergence of smooth linear statistics S/ 
for any determinantal process with correlation kernel . 

Definition 4.9. In the sequel, the quantity 2tB(^ / f{t)^dt will be called the Poisson component of the 

Jr 

cumulant and we will use the decomposition lim C^n [S/]=2rB((, / /(t)”dt + 0((, ^[/] where 

N^oo Jjj 


J •'»< i=l |m|=n 


X . 


(4.33) 
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This name is motivated by linear statistics of the Poisson point process whose cumulants are equal to t' 
where r' is the intensity. 


fitTdt 


Corollary 4.10. Consider the determinantal process with correlation kernel ^ which satisfies the con¬ 
ditions (IQl-TTKll and let f G The random variable S/ converges in distribution as N ^ oo to a 

random variable where t' = jdr and whose cumulants are given by 




(4.34) 

Proof. We can estimate the growth of ©((, ^.ilf], (I4.33|) . and the Poisson component separately. We start by 
giving an upper-bound for the constant . Recall equation (14.2911 . 


P n 

= n! / E M(m)E„ d 

.'-1 


X . 


I m| —n 

Obviously for any x G M", E„ [a;cr(s^,)] < Xn and if we use formula (I4.32|) . 

1 


it nit 1 /* 

TT $(a:i)E„ [xs ] d"a; < / TT ^{xi)xnd'^x = - -— / (1 - '^{x)Y~^^{x)xdx . 




Moreover, since 0 < $ = —and 0 < ^ < 1, we have for any n > 1, 

p poo poo 

/ (1 — 4'(a;))"“^$(a:)xdx < / ^{x)xdx = / ^{x)dx . 

Jr Jo Jo 

On the other hand, if we use that xi < E„ [a;cr(s^)] and apply the same method, we can show that 

[ TT $(xi)E„ [a;sj d^-x > - — [ 'i/{x)'^~^^{x)xdx > . [ (l - '^{x))dx . 

Jr’^ (,n—ijljR (^n—ij! 

These estimates show that there exists a positive constant C which only depends on the shape 4' such that 
for any 1 = 1,..., n, 




i=l 


< Cn . 


Using the estimate of lemma 33] and (14.291) this implies that | B(^ | < C{n 1)!2 
Hence, for any n > 2, the Poisson component is bounded by 


■n— 1 


b:^ / fitrdt 


<C"(n+l)!(2||/||oo)" 'll/IUi . 


(4.35) 


In the second half of the proof, we estimate the growth of ©^ [/]. Applying the upper-bound (14.231) in 

formula 1|4.33|) . we see that 


©$,r'[/]| < 2n! E |M(m)| f /(n*) (|mi| H - h|u„|)(i” f ^<^{xi)d"-x . (4.36) 

|m|=n ■’^0 •'*< i=l 

By symmetry f TT ^{xi)d^x = ^ ( f 4 ’(a;)da; ) = ^ and lemmas 33 ] and 331 provide bounds for the 

tJi 2 ri\ 

+ \\f\\mr ■ (4.37) 


other factors of the r.h.s. of (14.361) . We obtain 

|©^...[/]| < (n + l)!4^ 
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The estimates (I4.35|l and (14.3711 show that the limits of theorem 14.41 satisfy the criterion of lemma 14.81 for 
any choice of parameters r' > 0, dt G 5^, and / G i7o(K)- Hence they corresponds to the cumulants of some 
random variable which we denote by and 5/ ^ ’E.qf^r'f as TV —>■ oo. □ 


For the critical modified CUEs, according to formulae (I4.3H 1 and (14.41) . the parameter /3 = 1. Hence, 
theorem O follows directly from proposition 12.41 and corollary 14.101 Likewise, in the GUE setting, by 
formulae (ICTH. the parameter f3 = 1/4. Provided that 1/3 < a < 1, by proposition 11.91 we conclude that 
at the critical scale, a modified GUE with shape G 5^' converges in distribution to the random field t-/ 4 . 
In order to deal with all mesoscopic scales, we can use the asymptotic expansion of theorem 15.21 instead of 
theorem 15.31 Namely, if we combine formula (13.71) with the sine-kernel approximation (15.10|1 . we obtain for 
any scales 0 < a, 5 < 1, 


N-^Kl^{xN-\yN-^) 


1 

tN°‘ 


E 

|fe|<rAr“ 


$ 


/ k \ 

sin 

{N + k) (f 

f TT 

(2 ^N+k / 

-F 

(TTyN^/'^-^\ 
\2 ^N+k ) 

)] 


7r(x - y) 


-h O (N-^). 

N^oo 


(4.38) 

The r.h.s. of (14.381) is not a translation-invariant, so we cannot defined its Eourier transform. However, it is 
related to the kernel (11.21L by a change of variables and we can exploit this fact to compute the limits 
of critical linear statistics of the modified GUE at any scale, including the regime 0 < a < 1/3. 


Proposition 4.11. Let 4' G 5^', / G Co(K.), and 0 < a < 1. For any n> 2 


lim 

N—^OO 


[S/a] 


lim C 
N—^OO 


tN 


^9n] 


where 

ffjv(x) = /0iV“G(^)^ (4.39) 

and the function G is given by definition 15.11 

Proof. Let T = (logfV)^ and we assume that supp(/) C [—L,L]. 

Observe that for any |fc| < TiV" and any x, y G [—2L, 2L], a Taylor expansion give 


F X 


J^l/2-a 

\/ N + k 


-F\y 


’ y/WTk 


]ffrk (i^) - (iw)} U2„ ■ 


Thus, taking 5 = a in equation (14.38|) . we get for any 0 < a < 1, 

1 




tN° 


E 


/ fi \ sin 

^N{N + k) (F(l^)-F(f^)) 

) 

7r(x - y) 


+ O (N-°‘) 

N^oo 


\k\<rN° 

(4.40) 

where the error term is uniform for all x,y G [—L,L]. Eollowing the proof of lemma [221 this approximation 
implies that for any any composition m, 


Tr[/< 




t^N 

’a -^^,c 


E 

/cGZ^ 1=1 
IUI<rJV“ 


Jb. 

tJV^ 




Sin 


^N(N + k)(F( 


TT 

2 N° 


)-f (§¥")) 


7 t ( 3 -, -Xj+l) 


— (TX + 


O (iV-“), (4.41) 

N —¥oo 


where x^+i = Xi. There exists iV^ G N such that for all N > Ni^, we can make the change of variables 
yj = Tr~^N^F (f ;^) in the integral (14.411) . Since 0 < F'{x) < 2, this change of variables maps the interval 
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\xj\ < L to some subset of \yj\ < L for any N > N^. Hence, if we let gN{y) = f {^)) and r]{k) is 

given by (I4.3H H. we get 


/ YlfixjY 

'[-l.lv 


Sin 


VNiN + k) {F (1^) - F (f ^)) 




7T{Xj - Xj+l) 

n n;,(v ~ Sx 

LL9N[yj) -G(^)) 


J[-L.LV^^, 

A Taylor expansion gives for any y,z £ [—4L, 4L], 

G' {yN-^y^ iV“ {G{yN-y - G{zN-y} = {y-z){l + {{y - z)N-y} . 


(4.42) 


This implies that for any |fc| < rA^“ , 

G' (Yy) sin [2TTy{k) {y - z)] _ sin [2TTy{k) {y - z)] 

N^G {§L) - G {^)) - yy-z) 

Hence, if we combine formulae (14.411) , (14.421) and (14.431) , we have proved that 


+ o (N-y . 

N^oa ^ ' 


(4.43) 


i 

= E n 

feGZ^ i=l 
|fcil<rAf“ 

By (11.211) , we can write this equation as 

Tr[/r<---/r<] = Tr[5]^^L^,,...57iE] + „0 (fV-“) , 


and the proposition follows from formula (11.331) . □ 


_kj_ 

tN^ 


[-iT]%=i 


X{9N{.yy 


. sin [2Trri{k) {y — z)] 
Tr{y - z) 


O (N- 

N^oo '■ 


By proposition HTI] to complete the proof of theorem II. 61 it remains to extend the argument of theorem 14.41 
to deal with test functions of the form (14.391) . As we can see from the following lemma, such functions depend 
mildly on the density N and it is not difficult to finish the proof, see proposition 14. 131 below. 

Lemma 4.12. Let f £ G^(R) with support in [—L, L] and 0 < i5 < 1. For any N > ( 2 L)^F^ ^/jg function 
gN{x) = / has compact support in [—L,L]. Moreover, we can check that 

\\ 9 N-f\\oo= O {N-A and -/||^i = o(l) . 

N^ao 


Proof. By definition 15.11 0 < F' < 2 and the map x i—>■ ^N^G {j^) is a dilation. Therefore, when N > 
(2F)^/'^, the map gN is well-defined on [—L,L] and gN{x) = 0 for all x £ [—F, F]\supp(/). By continuity, 
we can assume that gnix) = 0 for all |a;| > L. Hence g^ £ Go(R.) with supp( 57 v) C supp(/). 

Then, by Lipschitz continuity of /, for any \x\ < L, 


IotW-/WI<cjv'|2g(^)-^ 


< C’L'^N 


-s 


where we used that G is smooth with G(0) = 0. This implies that 


||5Af-/lloo< [ \gN{x) - f{x)\dx = 0{N A- 
J[-L,L] 
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Similarly 


\9M - f{x)\ < ll/'lloo \2G' (- l| + |/' (^)) - f'{x) 

Since G'(0) = ^ and /' is continuous, lim \g'pf{x) — f'{x)\ = 0 for all x € [—L,L]. By (I1.29P and the 

N—^OO 

dominated convergence theorem, we conclude that 

\\9N - fWm = ^ ^ \9Nix) - f{x)\^dx 0 


as N 


n 


Proposition 4.13. Let / S Go(K), 0 < a < 1, and 'P G . If rj is given by (I4.3k zl and is given by 
(|4.39l) . then for any n>2, 

lim G2« [Esn] = C" [E^.r/4f] ■ 

Proof. We can take / = in formula (j4.8l) . In particular the estimate (14.181) is still valid for the test 
function g^- Let |m|i = |mi| + ■ ■ ■ + |Mn| and the function T^(u) be given by (14.141) . We want to prove that 


lim / 1; „i_oi 

v“|> 2 I 


and 


lim 

N—^OO 


n ^u=o, 

i 


(4.44) 


(4.45) 


Like (I4.19|l . these limits imply that 


lim G^n [= 5 jv] = 2 lim f K i TT 5 Ar(M*) i T)(r(M) d” 

N^OO N^OO J]gn I J 

= 2 lim f K < TT f{ui) \ T^iu) d^~^u , 

and the rest of the proof is identical to the proof of theorem 14.41 To complete our argument, it remains to 
show (14.441) and (14.451) . First observe that 

^ [ ^^hu^\yN]^\\^hT[\9N{ui)\d'^-\ . (4.46) 

JR" I 2 / ^ yjjn \ 2n f ^ 

Let An = |u G M : |u| > | and define the function qn by its Fourier transform qn = ^An9n- Then by 

(14.461) and lemma H751 we have 

l{|„|>i^}l“li n \9N{ui)\d^-^u < 71^2"-^ (IldArlloo + lldwllffi) (||5 a/||oo + IldArllffO"”^ ■ (4-47) 


(4.48) 


n—1 


A change of variables yields 


Hence 


\f{z)\dz. 

ffArIloo < II/IIlg and by definition || 5 jv||hi = ||/||hi, so that by (I4.47F 

/ l“l n \9nM\ d-\ < (IlgvrIU + \\qN\\m) (||/IUi + \\f\\m) 

■1 Kn I 2 / 
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Obviously ||giv||oo —>■ 0 and to conclude that (14.4411 holds it remains to estimate ||gAr||_f/i. The main observation 
is that 

r . 2 


qN{v) - 1 . 4 n /(^) \wd'>^ < \\9N - fWm ■ 


Then, by the triangle inequality, 


lkiv||l/i < 2 


llffw - /II 


2 





The first term converges to 0 by lemma 14.121 and so does the second term by the dominated convergence 
theorem since / G H^(R). We conclude that ||(Zjv||hi —^ 0 and we have computed the limit (14.441) . 

The proof of (14.451) is very similar. We observe that 


]^g7v('Ui) - n/(u*) = (5Af(%) - /K)) • 

i i j=l i>j i<j 

By (14.241) . there exists some constant C which only depends on n and 4' such that |T^(u)| < Cjl + |u|} 
and we obtain 




n « 

|T^(u)| (T-^u <CY 9N{uj) - f{uj) 


j = l •' *0 




i>j 


i<j 


{l + |M|}d" . 


Thus, by lemma IT751 and the facts 


< ll/llii and bwllffi = WfWm, we get 


/ 




Ui, 


T^(u)| < Cn^2"-^ ||g^ - /|U + llffiv - f\\m iWfh^ + \\f\\m) 


\ n—1 


(4.49) 

Lemma l4.12l implies that the r.h.s. of equation (I4.49P converges to 0 as —>■ oo and the limit (14.4511 follows. □ 


By proposition l4.11l and 14.131 if / G Co(M), for any 0 < a < 1 and any n > 2, 

lim [EU] = C" [E^,r/ 4 f] ■ (4.50) 

N—yoo -''-'I',a 

In the proof of corollarv l4.10l we have proved that the sequence (C" [S,j, ,-/ 4 /]) given by theorem [4.4l satisfies 
the condition of lemma IT51 This implies that, considering the determinantal process with correlation kernel 
the random variable 

S/„ ^ (4.51) 

as iV — >• oo. In order to complete the proof of theorem ll.61 we use a density argument to extend (14.511) to all 
test functions in 


Proof of theorem 1 1. dl First observe that for any z, y G K, 

I <4|a._y|2 ^ 


By Chebychev’s inequality, this implies that, if X and Y are random variables defined on the same probability 
space, for any ^ G K, 

|E I < 4|C|v'Var[X-y] . 

For critical linear statistics of the modified GUE, using the estimate ()3.16|) which is valid when S = a, this 
implies that there exists a constant C > 1 V 'J^^- such that any test functions /, G 7Jg^^(K.), 


I < C\f\^\\f-h\\l, + \\f-h\\l,,, . 


(4.52) 
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Moreover, formula (jA.4ll implies that under the same assumptions, 


|E I < C\^\^\\f-h\\l, + \\f-h\\ 


2 

// 1/2 


(4.53) 


By the triangle inequality. 


1 




< 


_)_ |]g _ gi?S;[,_.r/4/j | _j_ 



If we suppose that h S C'o(K), by (I4.51L the last term in the r.h.s. of (14.5411 converges to 0 as iV —>• oo. Thus, 


using the upper-bound (j4.52p and (|4.53p . for any / S and C S M, 


lim 

N—^OO 


(4.55) 


Since, the space Cg is dense in i?g^^ with respect to the norm ^ II-Hi.+ 11-11?,./.. the r.h.s. of (14.551) is 
arbitrary small by choosing h € Cq (K) appropriately, and we conclude that S/q, T/ 4 f as —>■ oo. □ 


4.3 Properties of the random process 

In this section, we study the random variables which arise from the limit of linear statistics of the 

critical modified ensembles. Because of the complicated structure of the cumulants in theorem ll.61 we cannot 
get much information about the random fields except that they are not Gaussian. However, as we expect 
from figure [U we recover Gaussian fluctuations in both limits r —>■ oo or r = 0; see proposition 14 .1 5l We can 
go a bit further and compute the Laplace transform of the Poisson component; see formula (14.661) below. For 
the field corresponding to the MNS shape function tjj{x) = (l-|-e““)“^, it turns out that this component 
is Gaussian, see equation (14.691) . This special property should be related to the fact that the MNS model 
arise from a Grand Canonical ensemble. A natural but difficult question is whether the sequence given 
by (14.331) also corresponds to the cumulants of some random variables, so that the field would be the 
superposition of white noise and an independent non-Gaussian process. 

Proposition 4.14. If the shape 4 € S' satisfies the condition B(^ ^ 0 for some n > 2. Then, for any t > 0, 
the random process of eorollarv \A.lQ\ is not Gaussian. 

Proof. It is clear from the definition p.38|) that lim G™(u, a;) = 0 and it follows from (14.331) that for any 

T—^OO 

n>2 and / G i7g(K), lim (/) = 0. Hence, by formula (14.341) 

T—^OO ' 

C" [E4,.rf] = 2rB^ /" fifTdt + o(l) . (4.56) 

jR r^oo 

Thus, the Poisson component dominates at large temperature and the random field of S^i.r/ is not Gaussian 
since there are test functions such that C"[Sif.T/] 0 whenever B((, ^ 0. This observation is actually valid 
at any temperature t > 0 because of the scaling property of the cumulants. By definition, G™(m, a;) = 
rG™(u/T, a:) and the change of variables Ui = rvi leads to 


© 


n 


[/] = -2 




d^-\d^x . 


Hence, by (14.341) . the random variables ^^,t/ and S®,i/(y) have the same distribution. □ 

Proposition 4.15. For any function f € Hg(]R), the rescaled random variable ^^f converges in 

distribution as t ^ oo to a Gaussian random variable with variance B|, ||/||^ 2 - On the other hand, 
converges in distribution as t ^ 0 to a Gaussian random variable with variance ||/||^i/ 2 - 
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Proof. When r —>■ oo, the asymptotic of the cumulants of the random variables are given by equa¬ 

tion (14.561) and 

C” =t1-"/2b^ [ fitrdt+ o (r-"/") . 

Jr 

H6I1C6 

lim C" [S^,.(r-i/V)] =<52(n)B|||/||i. . 

T—^OC 

Taking the limit as r —?> 0 is more sophisticated. We shall see that we recover the cumulants of the sine 
process given by lemma 2 in [35) and the cancelation follows from the Main combinatorial Lemma. We fix 
some composition m of n and some vector u G Mg ^'^el we will look at the symmetries of the function G™. 
By definition (11.3811 . 

G™(u,a;)= max|A'J(u)| (4.57) 

crGS(n) 

where s = (cr), (11.371) . The important fact is that this expression becomes independent of the variable x. 

In the sequel, we will denote G™(m) instead of G^{u,x) and we define 

Z\ (r^) (^1 ‘ ‘ ^rni j “t“ * * ' “t“ Ujn 2 7 ■ * ■ ; T ‘ ‘ ' “t” , 0) • (4.58) 

By definition (11.361) . for any u G Kg, we have 


{ Urns + l-\ - i-Urrii if S < i 

Unis + l + • • • + Wn + Ml + • • • + Umt if i < S 
0 if * = s 


For any s = 1,... ,i(ni) we let Tr^ G S(n) be the cyclic permutation given by 

Trs(i) = Wis -b i mod n . 

Then, we see that {A™(u) : i = 1, •••,£} = {A™(7rsu) : i = 1, • • • ,£} and, by (j4.57l) . we obtain 


G™(m) = Y ™ax{^r(7rsM)} , 

fTGS(n) “ 

Y G™(7ru)=n! Y . (4.59) 

7rGS(n) 7T£S{n) 

By dominated convergence, we can pass to the limit r —>• 0 in formulae (14.331) and (14.341) . By (I4.57L the two 
integrals decouple and we get 


lim C” 

t \,0 


[E7,yf]=-2 


I m| —n 


(4.60) 


This limit is independent of the shape 41 and it will be denoted by C" [^o/] ■ If we use the notation (14.581) . 
Soshnikov’s Main combinatorial Lemma reads for any u G Kg, 

E E M(m) max{A™(7ru)} = —(52 (m)|ui| . (4.61) 

7rGS(n) |m|=n 

Next we symmetrize formula (14.601) over all permutations of rt, by equations (14.591) and (I4.6ip . we conclude 
that 

C” [Sg/] = 62 in) f f{u)f{-u)\u\du . 

Jr 
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This shows that the random field .no is Gaussian with covariance structure (/, g)ffi/ 2 . 


□ 


In the proof of proposition 14. 141 we have seen that r is just a scaling parameter. Therefore, in the sequel, we 
will assume that r = 1 and write G™ = G™, etc. 


By definition 14.91 the behavior of the Poisson component of the field is encoded by the constants B^, 
(|1.391) . In the remainder of this section, we will compute the generating function of the sequence B^ and 
prove proposition 11.71 We start by a combinatorial lemma and a definition. 

Proposition 4.16. For any z,w€C such that — 1)1 < |l + w|. 


oo n— 1 






n\ 


1 + 


(4.62) 


Proof. By equation (11.4011 . 


1=1 


n — I 
k + 1 — I 


E 


ni ,...,71^ >1 

niH- \-ni—n 


nil - ■-nil 


So that if we exchange the order of summation between k and I, 


JFL 2z} yo^+'^Z^ — — 


n=lfc=0 


n\ 


, E 

n—1 1—1 ni,...,ni>l 

niH- ^ni—n 

oo n 

= E^’^E E 

n—1 1—1 ni,...,ni>l 

niH- ^ni—n 


(-l)m 

n—1 

V ( 

' n — l 

ni! • • • ml 

^ 1 

^k + 1 — 1 


k—l — 1 


(-l)m 

1 ^.1 

:w\l + 

w)^-‘ . 


W 


fe+1 


nil‘“ni\ 


Then, since E«“ E 

n—1 ni,...,ni'>l 

n\-\ - \-ni—n 

between I and n, we obtain 


1 


ni! • • • ml 


= (e“ — 1)^ for any a G C, if we exchange the order of summation 


EE^;j^ = E(-i) 

n—11—1 




1 + w 




□ 


This proves (14.621) using the identity —7’ 1^1 

In particular, we can deduce from proposition 14 .1 6l that the triangular array 6^ has the following symmetry, 
for any for fc = 0,..., n — 1, 

. (4.63) 

Definition 4.17. A shape function G is called symmetric if its derivative —$ is even. In other 
words, if it satisfies for all x G M, 

1 — 4'(—x) = 4'(x) . 

Then, for any symmetric shape, the map 


n—1 


X !->■ $(x) Ew^)'=(i-'i^(^)) 


n—l — k 
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is even when the index n is odd. Thus, for any m > 1, 


T:)2m+1 


P n—1 

k=0 


(4.64) 


If we substitute w = and z = ^(1 — 'I'(x)), for some |^| < e into (|4.62|1 . we get 


£ S E - fMr'-' - 

n=l 


1 + 4'(x)(e5 — l) 

Integrating both sides, by definition (11.391) . we get the generating function 

— 1 


E = / x4.(x) 

n=l 


1 + 'I'(a:)(e? — l) 


dx . 


Hence, the Laplace transform of the Poisson component of the random process is given by 

e«/(‘) _ 1 


^ pn r pp 

/ f(trdt= // x$(x) 
n- JR JJ 

RxB 


1 + '^{x){^e^^d) — l) 


dxdt 


(4.65) 


(4.66) 


for any test function / G C'o(R). We end this section by the proof of proposition 11.71 We can interpret 
equation (|4.69|1 below as the fact that the Poisson component of the field is a White noise. 

Proof of proposition [13 Letting 4'(a;) = tp{x) = (1 + e“) ^ in equation (|4.65p and p = — 1, we get 


E|fb; = 


= P 


f 


(l + e^)(l + e^+p) 
“ logt 


dx 


0 (l+t)(l + t + p) 


dt . 


We want to prove that for any p > 0, 


i 


0 (l + t)(l+t + p) 2 


Indeed, since p = — 1, formulae (I4.67|) and (j4.68|) imply that 


°° era 1 e2 

E|Bj = -iogV)4, 

n—1 


(4.67) 


(4.68) 


(4.69) 


for any |^| < e To prove (I4.68|) . we can differentiate both sides with respect to the parameter p and we 
see that it is enough to prove that for any p > 0, 


f 

JO 


lo (1 + i + p)^ 

An integration by part yields for any e > 0, 


logt ^^^log(l + p) 


1 + p 


(4.70) 


i: 


\ogt 


rdt = 


—eloge 


(l + t + p)2 (i + p)(i + e + p) 

Taking e —>■ 0 in this equation gives (14.7011 and then (I4.68|) follows. 


log(l + e + p) 
1 + p 


□ 


For general shape function, we do not know whether the Laplace transform of the Poisson component (14.661) 
arises from some random processes. 
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4.4 The third and fourth cumulants 


Proposition 11.71 implies that we cannot apply proposition 14.141 to rule out that the process is not Gaus¬ 
sian. The goal of this section is to prove that indeed it is not Gaussian. Our first attempt is to compute 
the third cumulant of the random variable S^/, but it turns out that it vanishes for any test function; see 
proposition 14.191 Consequently, we construct a test function y G 5(M) such that ^ 0. 

In the r.h.s. of (14.331) we can symmetrize the functions G™ with respect to permutations of the variables 
Ui in order to simplify the expression of the cumulants of the random variables S^./. There are even more 
simplifications possible using e.g. the constraints ui + ■ ■ ■ + Un = 0 and the DHK formulae; see remark 14.201 
below. In general, every step is elementary but given the complexity of the functions G™, (ll.38p . it becomes 
really difficult to systematically simplify equation (14.3311 . Nevertheless, to provide our example, we have to 
proceed for the 3'^'^ and 4*^ cumulants. All these computations are performed in appendix |B] For the 3'^'^ 
cumulant, we obtain the following formula. 

Lemma 4.18. Let ^ be a symmetric shape (definition \A.17\ and define the function ro : —>■ by 

'^(vi,V2) = ['Cl]'*’ + ['C 2 ]''' + + V2\^ — 2max{0, v\,vi + V2} ■ 

For any function f G i7o(R); 

C^ [S^/] = 12 f f ‘f’(s) ff A?(s + zi)A>(s + zi+Z 2 )'!a 7 (ui-Z 1,U2-Z2) d^udsd'^z , (4.71) 

I o' I IR 

° * '' ( 0 . 00)2 


where dfu = duidu 2 and it is understood that U 3 = —ui — U 2 - 

Proof. Appendix [B) □ 


The important feature of formula (14.711) is that the functions / and $ are coupled by a function vu which 
only depends on the variables Ui — Zi. Moreover, it follows from the DHK formula (14.7211 that this function is 
anti-symmetric and we will show that consequently the 3'''^ cumulant vanishes for any test function. We shall 
see that the higher-order cumulants do not enjoy such symmetries. This means that for the critical modified 
ensembles, there is no cancellation such as Soshnikov’s Main Combinatorial Lemma. 

Proposition 4.19. Assume that the shape At G ^ is symmetric. For any function f G i7g(R), we have 

[E^f] = 0 . 


Proof. For any ru G R^, we have 


max{0, vi,vi + V 2 } + max{0, V 2 , vi + V 2 } = [ui]+ -I- [^ 2 ]^ + [ui -I- ^ 2 ]^ 


(4.72) 


It follows that another expression for the function w is given by 


w{vi,V 2 ) = max{0, V 2 , vi + V 2 } — max{0, vi,vi + V 2 } 
In particular zu(vi,V 2 ) = —■a7(v2,vi) and it follows that 
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<i>(a;)<i>(a: + Zi)^{x + zi + Z2)w(ui — zi,U2 — Z2) dxd^z ( 4 . 73 ) 

Rx (0,C5 o)2 

= JJJ ^{-y - zi - Z2)^{-y - Z2)^{-y)w{ui - zi,U2 - Z2) dyd'^z 

Rx (0,oo)2 

= JJJ ^{y)‘^{y + zi)^{y + zi + Z2)w{ui - Z2,U2 - zi) dyd'^z 

Rx (0,cso)2 

^{y)^{y + Zi)^{y + 2:1+ Z2)w{u2 - Zi, Ml - Z2) dyd'^z . ( 4 . 74 ) 

Rx (0,oo)^ 

At first, we made the change of variable y = —x — zi — Z2- In the second equality we used the assumption 
that 4> is symmetric, and in the last equality we used the anti-symmetry of w. Equation (14.741) shows that 
the integral (I4.73P changes sign under permutation of the variables ui and U2- Because of this fact and the 
symmetry of formula (14.711) . the 3'''^ cumulant of the random variable vanishes. □ 




Remark 4.20. The Dyson, Hunt, Kac (DHK) formulae are the following remarkable identities. For any 
n> 2, 

"1 + 

^ ^ max “t“ 5 ■ • ■ 5 ^7r(l) 4“ * * * 4“ 1) : b} ^ ^ ^ ^ y [^7r(l) 4“ * * * 4“ ^7r(/)] ■ (^■^^) 

7rGS(n) 7rGS(n) 1—1 

When n = 2, this gives equation (j4.72D . We refer to Simon’s book section 6.5 for a proof of (14.751) and 
an application to the Strong Szego theorem. Actually, it is also possible to give a proof of Sosnhikov’s Main 
Combinatorial Lemma (I4.61|) based only on these formulae. We can also apply (14.751) to the cumulants of 
linear statistics of the modified ensembles but this only leads to partial simplifications. 

To compute the 4*^ cumulant, we make the change of variables x G i—>■ (s, z) G R x Ri). given by 


S = Xi Zi = X2 — Xi Z2 = X3 — X2 Z3 = Xi — X3 


in equation (14.331) . We get 


«>/-(/) =-2 


[ e{z) [?ft\Y[f{u.) 

^4 1-1 


|m|=4 


(u, z)d^ud^z 


(4.76) 


where G™(m, z) is the image of a; !->■ G’"^(u,x) (observe that it does not depend on the variable s) and 

0(z) = / 0(s)(()(s 4-Zi)((i(s 4-Zi 4-Z 2 )(^(s 4-Zi 4-Z 2 4-Z 3 )ds . (4.77) 

It is worth noting that, since the function (p is even, we have 0(z3,Z2,zi) = 0(zi,Z2,Z3) but no further 
symmetry. Unlike lemma 14.181 we cannot obtain a compact formula for the 4*^ cumulant because there is not 
enough symmetry. Thus, it is simpler to compute the value of the functions G {u, z) at some given points u 
and deduce that ©^(y) 0 for some test function y which is sufficiently concentrated around these points. 

The technical result that we need is given in the following lemma. 
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Lemma 4.21. For any z € up to the permutation of zi and Z 3 , we have 

Y, 1 : M(m)G"(„, ^) = 24 (4[1 - 4,1+ + 2[1 - 4. - 4,]+ + ^ 

JJIH-1-«4=0 |m|=4 

«ie{-i,i} 

+ [2 — zi — Z 2 ]'^ ~ [2 — zi]^ — 2 max{0, 1 — zi, 2 — zi — Z 2 } — 2 max{0, 1 — zi, 2 — zi — Z 2 — za}^ . 
Proof. Appendix [BJ □ 


It follows from lemma H. 2 II that 


H / 0(z) ^ M(m)G"'(z;,z) d3z = 0.29... . (4.78) 

„i+...+4,4=0-^(0.°°)^ |m|=4 

fiGl-l,!} 


The integral can be performed analytically or numerically using Mathematica. To complete our argument, 
we also need the following approximation lemma. Its proof is rather straightforward and for completeness it 
will be given after our example. 


Lemma 4.22. Let 


lim 

£->■0 


g{x) = e and'Qi{u)= / 0(z) M(m)G™(M, z)d^z. For any v , we have 

|4„|=4 

j ^ g ...g 2t(u) d^u = (5o(ui +- h U 4 ) . 


We let y{x) = 2e cos(27ra;) for some e > 0. Then 

and by lemma 11.221 we have 

[[ Y[Mui)e{z) M{m.)G’^{v,z)d^zd^u=^ ^ 2t(u) . (4.79) 

'^R^(0.oo)3 * IH=4 

The r.h.s. of equation (I4.79P is given by (14.781) and by equation (14.761) 

lim ©+(?/) = —0.29... . 

e —>-0 

Since the constant = 0 by proposition 11.71 we conclude that, if the parameter e is sufficiently small, 
G^ [^^y] ^ 0 and the linear statistics E^y is not Gaussian. 


Proof of lemma \J^^ Let us fix u G and let r{v) = ui H-+ U 4 . It is easy to see that the functions u !->■ 

G™(m, z) are Lipchitz continuous with respect to |u|oo with some constant which can be chosen independently 
of z G Then, the function 21 (m) is also Lipschitz continuous on R^. A change of variables yields 


^-4 


Ul - Vi 


9 


9 


U 4 — V 4 


2 l(u) d^u = 


g{wi) ■ ■ ■ g{w4)i2i(v + ew) d^w 






















If we let € M{e ^r(v), ^), it is easy to see that 


So that if + • • • + t’4 = 0, i.e. r(v) = 0, then E 

_4 / /ui - vi 


g(wi) ■ ■ ■ g(w4) (fw = -^E 




(4.80) 




= for any e > 0 and it follows that 


lim e 

e-)-0 


/’ 


■9 


U4 — V4 


2l(w) d'^u = -2t(r;) 


On the other hand if r{v) ^ 0, by equation (14.SOL 

lim / g{wi) ■ ■ ■ g(w 4 ) d'^w = Q 


r{v)/e 


□ 


5 Asymptotic formulae 


In this section, we give some background on the asymptotics of the Hermite functions and the GUE kernel. 
Asymptotics of orthogonal polynomials have been extensively studied, see for instance m or m- For appli¬ 
cations to statistics of Orthogonal Polynomial ensembles, one is generally interested in uniform asymptotics 
of the Christoffel-Darboux kernels. Based on the non-linear steepest descent method m, for a large class 
of potentials, sine-kernel asymptotics have been established in [12] at the microscopic scale and recently 
extended to mesoscopic scales in [25]. We review some aspects of the theory for the Hermite functions: 


K{x)=e^ ^ 


y^n!2"0r cJa:" 

Definition 5.1. We define on [—1,1] the functions g{t) = 2\/l — and 


(5.1) 


rx 

F{x) = / g{t)dt = arcsin(a:) -I- x\/l — . 

Jo 


The map F is a dijfeomorphism from |a;| < 1 to |a;| < ^ and we let G be its inverse. 

The Hermite functions have the following bulk asymptotics, for any 0 < 7 < 2/3 and for all |x| < 1 — 

1/4 


hn{V^x) = 


n(l — x^) 


1 


I cos ( — nE(x) — i arcsin(x) ) -I- O (n ^) I , (5.2) 

( y 2 2 J n—)-oo ' ' j 


where the error term is uniform. While the Hermite functions have oscillatory behavior inside the bulk, they 
have exponential decay outside. For all 0 < 7 < 2/3 and all |x| > 1 -I- 


hn{V^x) = 


\/v^SnTT Y^l — (|a;| -I- y/Ufi — 1) 4 




(5.3) 


where FI is an even function which is defined for all x > 1, by 

/ X 

g(t)dt = |x| \/x^ — 1 — log(|x| -I- \JX^ — . 

Note that in the Hermite case, using an integral representation, formulae (15.21) and (15.3|) can be obtained by 
the classical steepest descent method. At the edge, a precise asymptotics is also known but, since we are 
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only interested in bulk linear statistics in this paper, we won’t need any precise estimates and we will use 
instead the uniform bound 

||/ln||oo < CH , (5.4) 

where ch is a universal constant and the exponent is sharp. Moreover, if restricted to a fixed interval in 
the bulk, the error term in formula (15.21) can be improved. Namely, for any 0 < A < 1, uniformly over all 
X £ [—A, A], we have 


hn{-\/^x) = 


1/4 


1 


^n(l — x^) 

For the quadratic potential, we define the Christoffel-Darboux kernel 

N-l 


Icosfn - nFix) -arcsin(a;) ) + O (n 1 . (5.5) 

( \ 2 2 J n-loo J 


kU-,v) = i:' h,{^)K(y) = 

^0 V 2 x-y 


and the Wigner semicircular law, for any |t| < \/2, 


TT 7r\/2 v2 


(5.6) 


(5.7) 


At the microscopic scale, it is well-known that we get the sine kernel in the limit, 


lim ATcd 

y/m ^ V 




sin[7rgse(a^o)(^ - Q] 

- C) 


see e.g. [12]. This asymptotic has been extended to all mesoscopic scales in [25]. In [26], we give an elementary 
proof of theorem 15.21 which is based on applying the classical steepest descent method in the Mehler formula 
for the GUE kernel. In the sequel, L is some arbitrary large positive constant and |a;o| < V^- 


Theorem 5.2. For any —1/2< A < 1/2, we have the asymptotic formula 

sin A 


(yNxo + iN’', VNxo + CN^) = 

uniformly over all C S [—L,L]. 




Ke - C) 


-F o 

AT-j-oo V / 


This approximation takes into account the density of the Wigner semicircular law, i.e. the fact that the GUE 
eigenvalues are not uniformly distributed at the global scale. By definition 15.II and theorem l5.2l can be 
rephrased as 


sin 

A^A^d (yiVxo + VNxo + CN^) = — 


AI/2+a + tN-^/^+^)dt 


- C) 


-b o . 

N^oo \ ) 


(5.8) 

Note that unlike the sine-kernel, the kernel (15.81) is not translation-invariant. This gives raise to some com¬ 
plications to compute the cumulants of large scale linear statistics of the modified GUE, cf. proposition 
However, at sufficiently small scales, we recover the sine-kernel as a special case of theorem 15.21 

Theorem 5.3. For any —1/2 < A < 0, we have the asymptotic formula 

sin {N^^'^~'~^Trgsc{xo)(i - 0) 


A^A^d (VNxo + ^A^, VNxo + 
uniformly over all C S [—L, L], 




-b O (A^^) 
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Proof. By definition 15.11 


F{x) - F{y) = g 


x + y 


{x-y) + 0{\x-yf) 


If we let a;AT = ^ and yN = ^ Taylor expansion gives 


XN + yN 




and it follows that 


F{xn) - F{yN) = (C - C)N-^/^+^ + 0(1^ - C|iV-i+2^) 


Hence, by (15.71) . we have proved that 
sin 




V2 


V2 


(5.9) 


= Sin 


N^/^+\gs,{xom - O] + O (1C - . 

Af—>00 


When A < 0, the error term is converging to 0 for any x,y £ [—L,L], and if we plug this approximation in 
the formula of theorem 15.21 we obtain the asymptotics of theorem 15.31 □ 

Remark 5.4. In the special case xq = 0, the error term in dEll) is of order N 1+2-^ and the sine kernel 
approximation of theorem \5.3\ is valid in the whole range —1/2 < A < 1/6. 

The Christoffel-Darboux kernel is the same, up to a scaling, as the GUE kernel defined by p.4ll in 

the introduction. Namely, if we let A = | — i5 for <5 £ (0,1], by (15.61) . we can rewrite 


N-^K^{xN-\yN 


T-S 


) = 


-^xN^ 

72 





Thus, for any e > 0, if Xq = 0 and M = N F k iov some |fc| < by theorem 15.21 for all x, y £ [—L, L], 

we have 


sm 


N-^K^{xN-\yN-^) = 


M 


(^( 


2 VMTV■5-1/2 




2 


7r(x - y) 


+ O (N-^) 

N^oo 


(5.10) 


This formula holds at any mesoscopic scales. On the other hand, if we assume that 1/3 < 5 < 1 (cf. 
remark. since ysc(O) = by theorem 15.21 we get the following asymptotics 


N-^K^{xN-\yN-^) 


sin [rriV^/^ ^^M/N{x — y)] 
7r(x - y) 


+ O . 

N^oo ^ ' 


(5.11) 
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A Proof of theorem 11.5 


We give two different proofs of theorem If .51 First, we can use the ideas of section [XT] to compute the limit of 
the reproducing variance Vo{fs) at the critical scale 6 = a, see proposition lA. H and (IA.2I) . Since this result 
is not used in the rest of the paper, we will only sketch the argument. Second, we use (14.501) and compute 
r/if] using formula (11.1811 . Subsequently, we check that the different formulae for the critical variance 
are consistent and we apply them to the MNS ensemble. 

1 /2 

Proposition A.l. For any shape 'h G , any test function f G Hq' (R) and any 0 < a < 1, we have 




f{x) - f{y) 

2 

^(ix-y)\ 

x-y 


1 

H 


dxdy . 


Proof. By formula dSU) combined to the estimate (13.1011 and the approximation (14.401) . we see that 


Voifo.) = 


1 




E 


|fc|<r7v° 

b|<r7v° 


tN° 


$ 


tN'^ 


[-L,L]2 


fix) - fiy) 


x-y 


Sin 


FFTTF) (f (f - F (f i))] si„ [vWTTT(f (f - F (f i))] dxd, + , 


where the support of the test function / is included in [—■§, -j] and F = (log N)^. Note that there is another 
error-term of order N~°‘ coming from (14.401) that we neglected. Applying a trigonometric identity and the 
Riemann-Lebesgue lemma like in the proof of lemma 13.21 we see that 



[-L,L]2 


fix) - fjy) 
x-y 


2]^2o 


$ 


|fc|<rAr° 

|i|<rjv“ 




cos 


/ (fc - j)TTix 
2N^ 



dxdy + O {L ^) 

L—^oo 


as A —> oo. The sums converge to some Riemann integrals and by the dominated convergence theorem. 


lim 
N —^OO 


Voifa) 


J_ ff fix) - fjy) 
47r2 JJ x-y 

[-L,L]2 


<i)(M)<i)(u) cos 


7r(u — v)(x 

2^ 


y) 


dudvdxdy 


^ Cjf) 

- L 


(A.l) 


where the constant C{f) depends only on the test function / by (13.101) . To complete the proof, we can use 

2 


that for any gM., JJ $(u)4>(u) cos [27r(u — v)f] dudv = l>(^) and let L —^ oo in equation (lA.ip . □ 

As a consequence of lemma IXTl lemma IX2l and proposition lA.il we get for any 0 < a < 1 and / G iTQ^^(M), 


Jim^Var^.J/^] = ^ | f{x)^dx+^ JJ 


f{x) - fiy) 

2 


1 

H 


V 4r-i ) 


dxdy . 


(A.2) 


Remark A.2. Such a direct computation seems possible only for the variance of linear statistics thanks 
to the special structure of the reproducing variance Vq, en). A technical difficulty to compute the limit of 
the higher-order cumulants comes from the singularity of the correlation kernel ^ along the diagonal. 
Therefore, it is better to exploit instead the fact that the kernel yixen by (11.211) is translation-invariant 
and use Soshnikov’s method. 

We have seen in section ITX] that, up to a scaling, the modified ensembles have the same limit at the critical 
scale. Thus we have already obtained two expressions for the variance of the random variable 2,^ 
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formula (12.1211 in the Circular case, resp. formula (IA.2I1 in the Gaussian case. Theorem 11.61 provides a third 
expression, see (IA.4I) below, and we will check that all these expressions are consistent. By formula (11.181) . 


Var[S^..,/] = 2rB|||/||i.+2 / /(u) 

Jr 

By definition (11.361) . for any u G Rq, 


$(xi)$(x2)G^ {u,x) dudPx . 


(A.3) 


Xi <X2 


0 —U2 

U2 0 


A^\u)=(^ 

Then, by definition (ll.38p . for any xi < X 2 , 

G^ (m, x) = max{ 0 , U 2 — t{x 2 — cri)} + max{—U 2 — t{x 2 — a:i), 0 } = [|u 2 | — t{x 2 — cci)] 
Hence, for any function / G i?g^^(IR), by equation (IA.3I1 . 

Var = 2 tB| 11 / 11^2 + 2 /" /(u) [[ ^{xi)^{x 2 )[\ui\-Tix 2 - xi)]~''dudxidx 2 ■ 

JR JJ 


(A.4) 


Xi <X2 


We first check that this formula matches with the r.h.s. of (12.121) . Recall that $ = —4'' and, for any m > 0, 
some integrations by parts yield the identity 


<l>(t) f + s)[m — s]'''dsdt = f (l —'l'(t)) 

; Jo Jr 


4'(t + m) + - u) 


— 4'(t) ) dt . (A.5) 


By (12.121) and the definition of B^, we get 


lim (7^jvt/a] = 2 T / f{u) / 4'(t)(l — 4>(t + MT ^))dtdv 

N^oo P Jb 


= 2t 


[ f{u) f (l-'i’it)) 

Jr Jr 


4>(t + MT ^)+4'(t — MT 


dtdu 


= 2tB 


II J llL2 


2 r 


f{u) f ^(t) f + s) [\u\ — Ts]'^ dsdtdu 

Jr Jo 


which is equal to (IA.4I) . To check that formula (IA.4I) also matches with (|A.2I) . we use an argument that is 
similar to the proof of the identity (11.281) for the variance of the sine process. 

Lemma A.3. For any function f G any w > 0, we have 


1 

47r^ 


fix) - fiy) 


x-y 


cos (flirix — y)w')dxdy = / /(m) [|u| — w] du . 


(A. 6 ) 


Proof. By Plancherel’s formula, for any z G M, 

[ \ fix) - fix + z)\^ dx = 4 

JR 

Then, by Fubini’s theorem, for any w > 0, 


[ l/Wl^si 

JR 


sin^(7ritz)du . 


1 

47r^ 


fix) - fiy) 


x-y 


1 


cos (27r(a; — y)w)dxdy = 


fix) - fix + z) 




sin^(7rMz) 


cos ifhrzw^dxdz 
cos [flTTzw^dz j du . (A.7) 
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Moreover, by Residue calculus, one can show that for any w > 0 and u G 

sin^(7rMz) 


P 


■ cos(2TTwz)dz = 7r^[|u| — . 


The lemma follows by combining equations (IA.7I) and (|A.8|) . 


(A.8) 

□ 


If we take w = r(s — t) and integrate equation (IA.6I) against $(t)<i)(s) on the simplex {t < s}, we see that 
the so-called reproducing variance, see formula (lA.ll) . satisfies 




fix) - fiy) 


x-y 


$(t)$(s) cos (27rT(a; — y)(s — t))dxdydsdt 


/(wi) 


tOs 


‘l>(t)<l>(s) [|miI — r(s — t)]''"dxidx 2 . 


If we add up the Poisson variance, this implies that the r.h.s. of (IA.2I) and (IA.4I) agree. 

As an example, let us see what these formulae look like for the MNS ensemble. The MNS shape is ipP) = 
(1 -|- and an elementary integration gives 


f ij)[t)il — Ip{t + u))dt = f 

Jr Jo 

Then, by equation (I2.12L we get 


1 


0 (l + s)(s + e“) 


ds = 


1 — e~ 


Var [E^,rf] = / fiu) 


tanh(^) 

We can deduce the dual of this formula from equation (IA.2I) . We have 


du . 


Pit) = 


P(u) = - 


2tt^u 


so that 


cosh[t/2]2 ’ ' sinhpTT^M 

Var[S^,^/] fix)^dx + jj \fix)-fiy)\ 


sinh \2-k'^t{x — y)] 


dxdy . 


B Proofs of lemmas 14.181 and 14.21 

To prove lemmas [4.181 and 14.211 the strategy is to exploit the symmetry of equation (j4.33l) to simplify as much 
as possible the cumulants of the random variable S^/. To this end, we will use the following convention. 
Given two functions, we write / = y if there exists a permutation a G §(3) such that fiu) = gicru) or if 
f{u) = gi—u) for any u € K.^. For any sequence of real numbers, we will also denote 

mtixjui,..., Un} = max{0, ui ,..., Un} ■ 

Unfortunately the combinatorial structure behind the cumulants of the modified ensembles seems to be rather 
complicated and consequently the following computations are quite technical. 


Proof of lemma \jAB\ For any symmetric shape function 'k, we have seen that B|, = 0; see equation (14.641) . 
By corollary 14.101 this implies that 


[E^f] = ©|[/] 


= \ E {3G"+^(u,; 

|m|=n 


)+3G^+^(7/,x)-4G'' iu,x)\d^-^ud^x (B.l) 
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since M(2 + 1) = —| and M(IL^) = 2 by definition (11.321) . So it remains to compute the kernels and 

to prove equation (14.711) . By definition (I1.36L 


A^+^(w) 



Then, by definition (ll.38|l . 

= [ui — X 2 + + [ui — X3 + + [—ui — X 2 + + [—ui — X 3 + + [ui — X 3 + X 2 ]~^ + [—mi — a;3 + X 2 ]^ 

= 2( [mi - X2+ xi]^ + [mi - a;3 + xi]^ + [m - a;3 + X2]^ ) ■ (B.2) 


One can check that each term corresponds to a permutation in S(3) in the following order 123,132,213, 312, 231, 321, 
and at the second step, we used the symmetry between u and —u. By a similar argument, the contribution 
of the composition = (1,1,1) is given by 

( 0 —U 2 —U 2 — U 3 \ 

U2 0 -U3 j 

M2 + ^3 U3 0 / 

and 

Gi 3 (m, x) = max{u2 — X2 + xi,U2 + U3 — xs + xi} + max{u2 — X3+ xi,U2 + U3 — X2 + xi} 

+ max{—U2 — X2 + xi, M3 — a;3 + xi} + max{—^2 — xs + xi, M3 —X2 + X1} 

+ max{—M2 — M3 — X2 + xi, —M3 — X3 + Xi} + mtix{—M2 — M3 — X3 + xi, —M3 — X2 + Xi} . 


If we use the condition mi +M 2 +M3 = 0 and use the symmetry under permutations of the Ui’s, we can rewrite 
the previous formula, 

G]i 3 (m, x) = 6 mtix {mi — X 2 + xi,Ml + M2 — X3 + xi} . (B-3) 

If we combine equations (IB.21) and (IB.31) , we get 

3 G^''"^(m, x) + 3 G^^^(m, x) — 4 G^ (m, x) = 12 ^(mi — X 2 + xi, M2 — X3 + X2) . (B.4) 

Then, if we make the change of variables xi = x, zi = X2 — xi and Z 2 = X 3 — X 2 in formula (|B.ip . equation 
dEl implies that 

G^ [S^/] = 12 / i /(wi) i f ff ^{x + zi)^{x + zi + Z2)'cu{ui - zi,U2 - Z2) cPudx(fz , 

'J I I « M 'j'j 

'' * '' (0.Oo)2 

where it is understood that M3 = —mi — M 2 in the first integral. □ 


Proof of lemma \l.21\ We fix 2; G lRi)_. We will proceed exactly as in the proof of lemma H. 181 except that we 
will not give all the details. We will denote ± = +1 or —1 and we let 

Cl = mnx{l — zi, 2 — zi — Z 2 } + m^x{l — zi, 2 — zi — Z 2 — Z 3 } + max{l — zi — Z 2 ,2 — zi — Z 2 — Z 3 } 

+ 2[2 - zi]+ + [2 - zi - Z2]+ , 

C 2 = 2 ^[1 - Zi]+ + [1 - Zi - Z 2 ]^ + [1 - Zi - Z 2 - Z 3 ]'^'^ , 

Cs = 2^max{l - Z 2,2 - Z 2 - Z 3 } + [2 - Z 2 ]''' + [1 - Z 2 ]'^ + [1 - Z 2 - , 

C 4 = 4[1 — Zi]"*" + 2[1 — zi — Z2]~'’ , 

Cs = 4[1 - Z2]+ + 2[1 - Z2 - Z3]+ • 
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We will compute the values of G {v, z) for all compositions m of 4 and all points v = (±, ±, ±, ±) such that 
^ Ui = 0. The computations are not difficult but there are many cases to check. 

By definition (11.361) . 




0 —U 3 — rt4 

U 3 + U 4 0 


0 —Ui — U2 

lii + U2 0 


(B.5) 


Then, by defintion (I1.38L we can check that 

— ~ T [2 ~ '^2]'*' + [2 — 23]"*" + [2 — — 2:2]'*' + [2 — Z2 — 2'3]'*" + [2 — zi — Z2 ~ . 

We used that Ui + U2 = 0 or 2 and that > 0 to check which terms are a priori non-zero. Moreover, for the 
same reasons, 

G^^^+ - +-) = G^^^(+ - -+) + G^^^- + +-) = G^+^- + -+) = 0 . 


If we use the symmetry of the function 0, see (I4.77|) . under the change of variable zi -O- Z3 and that 
M(2 -I- 2) = —3, we can conclude that 

M(2 + 2) ^ ^ -24[[2 - zi]+ + + [2 - zi - ^2]+ + ~ ~~ ] (B.6) 

v-iA -l-i'4=0 ^ 

in the sense that if we replace the l.h.s. of equation (IB.61) by its r.h.s. in formula (14.761) . it does not change 
of the value of the integral. 

Let us continue with the compositions 3-1-1 and 1-1-3. We have 





-A^+^('u) . 


This expression depends on a single variable (say ui) and collecting the non-zero terms yields 

G (±, ±, ±, ±) = G^+^(±, ±, ±, ±) 

= 2([1 - zi]+ -I- [1 - 2:2]+ + [1 - Z3]''" + [l- zi- Z 2 ]+ -I- [1 - Z 2 - 2:3]’'" + [1 - 2;i - 22 - 23]+) . 
We can again use the symmetry of formula (14.761) and, since M(3 -I-1) = —2, we get 


M(3 -b 1) ^2 G (w, 2) -b G ^^{v, 2) i - 48f2[l — zi]+ -b [1 — 22]“'" -b 2[1 — 21 — 22]''" -b [1 — 21 — 22 — 23]''‘ 

«iH- \-V4=0 

fiGl-l.!} 


= -24(C2 + 2[1 - 2i]+ + 2[1 - Z2]+ + 2[1 - 21 - 22]+) . (B.7) 


Gonsider now the composition l"'^ = l-bl-bl-bl. By definition (11.361) . 


A^ 


4 


/ 0 -U2 

U 2 0 

U 2 + U3 U3 

\ W 2 + M 3 + W 4 M 3 + M4 


-M2 - M3 
-M3 
0 

M4 


— M2 — M3 — M4 
-M3 - M4 
— M4 
0 


\ 

J 
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If we look at all permutations in S( 4 ) and use the symmetry under permutations of the Mi’s, we get 


{u,z)= 4 (m 4 x{Mi — Zi, U1+U2 — Z1—Z2, Ui + U2 + U3 — Zi — Z2 — Z3} 

+ max{Mi — Zi — Z2 — Z3, Ml + M2 — Zi — Z2, Ml + M2 + M3 — Zi} 

+ max{Mi — Zi — Z2, Ml + M2 — Zi, Ml + M2 + M3 — Zi — Z2 — Z3} 

+ max{Mi — Zi — Z2 — Z3, Ml + M2 — Zi, Ml + M2 + M3 — Zi — Z2} 

+ max{Mi — Zi — Z2, Ml + M2 — Zi — Z2 — Z3, Ml + M2 + M3 — Zi} 

+ max{Mi — Zi, Ml + M2 — Zi — Z2 — Z3, Ml + M2 + M3 — Zi — Z2}) . 

So we can assume that G^ is given by the r.h.s. of this expression, then it is straightforward to check that 
G’^'*(+ H-) = S^m’axll - zi, 2 - zi - Z2} + m’ax{l - zi, 2 - zi - Z2 - Z3} + [2 - zi]+^ , 

G^"(- + +-) = G^"(+ - -+) = 4 C 2 , 

G^U+-+-) = <4 , 

G^\- - ++) = G^\- + -+) = 0 . 

So that, since M(ll^) = —6, 

M(ll^) G^\v,z) =-2A^2C2 + C4 + 2 max{l — zi, 2 — zi — Z2} + 2 max{l — zi, 2 — zi — Z2 — Z3} + 2[2 — zi]+ 

jjiH-|-i’4=0 

(B.8) 

If we combine formulae l|B.6l) . (IB. 71 ) and (IB. 81 ) . 


E E M(m)G'”(M, z) ^-24-| 3^2 + C 4 + 2m^x{l — zi, 2 — zi — Z 2 } + 2inax{l — zi, 2 — zi — Z 2 — Z 3 } 

i>lH-|-i;4=0 |m|=4 ^ 


1+ _L on _ ..„!+ _i_ on _ _ ..„!+ _L qro _ ^.1+ _i_ I?__|_ [2 - zi - Z2]+ + 


+ 2[1 — zi] + 2[1 — Z 2 ] + 2[1 — zi — Z 2 ] + 3[2 — zi] + 


(B.9) 


Finally, we look at the composition 2 + 1 + 1 , 


A^+^+^(m) 



-M 3 

-U3 - M 4 \ 

i j 

( 0 

-Ml 

-Ml - M2 \ 

M 3 

0 

-M 4 

- 

Ml 

0 

-M2 

\ M 3 + M 4 

M 4 

0 

' \ 

1 Ml + M2 

M2 

0 / 


and if we follow the same procedure, we can prove that 

g"+'+\+ + —) + G"+i+^(- - ++) = 2Ci + 2C2 + Ca , 

G^+^+i(+ _ +_) + G^+^+\_ + _+) = G^+^+\+ - -+) + G^+"+\- + +-) 

= 2^2 + C4 + Cs 


On the other hand 


^1+2+1 



-M2 - M3 

-M2 - U3-U4 \ 

, / 

0 Ml + M2 

M2 \ 

M2 + M3 

0 

-M 4 

- 

— Ml — M2 

0 

-Ml 

\ M2 + M3 + M4 

M 4 

0 J 

' \ 

, -M2 

Ml 

0 / 
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It is not difficult to see that, up to conjugation by some permutation matrix, we have 
This implies that G = G because such conjugation only changes the order of the sum over S(4) in 
the definition (|1.38l) . Similarly, we can check that the matrix A^+^+^ is also conjugated to A^+^+^, so that 
they give the same contribution to the 4*** cumulant. 

Since M(2 + l + l) = 4, putting all terms together, we conclude that 

M(m)G (u, z) = 24 + 3C2 + ^ + C4 + Csl ■ 

viA - i-V4=0 |m|=4 ^ ^ 


Observe that using the symmetry between Zi and Z 3 and the DHK formula (14.721) . we can show that 

Cl +1 ^ 2 (^[1 - ^l]+ + [1 - Z2]+ + [1 - Zl - Z2]+ + [2 - Zi]+ + + [2 - Zi - Z2]+ + ~ ~ 

and we get 


E E M(ni)G (v, z) < — 24< 3 C 2 T C 4 T Cs T 2[1 — + 2[1 — • 22 ]'^ + 2[1 — zi — Z 2 ]~^ 

vi-\ - \-V4=0 |m|=4 ^ 

I'iSf-l,!} f(m)=3 

+ 2[2 — zi]~^ ffi [2 — ^ 2 ]'^ ffi 2[2 — zi — -Z2]~*" + [2 — zi — Z 2 — •2'3]~*"^ ■ (B.IO) 


Finally, if we combine formulae (IB.91) and (IB.101) . a lot terms cancel but not all of them and we are left with 

iJiH-l-i'4=0 |m|=4 ^ 

— 2 max{l — zi, 2 — zi — Z2} — 2 max{l — zi, 2 — zi — Z2 — Z3} ) . 


Finally, if we make the change of variable zi Z3, 


Cs 4[1 — Z2]''’ + 2[1 — zi — Z2]^ , 


and we have proved the formula of lemma 14.211 □ 


Acknowledgement: We thank Maurice Duits for helpful discussions regarding his related works on 
mesoscopic linear statistics. 
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